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XXIV. BAND VIERTES HEFT 1956 


Zur algebraischen Mafsynthese ebener Kurbelgetriebe 
II. Mitteilung 
Von K.-H. Sieker 


1. Einleitung. Die algebraische Behandlung getriebesynthetischer Aufgaben bietet gegen- 
tiber geometrischen Methoden in verschiedener Hinsicht Vorteile. Die Darstellung in der Gauf- 
schen Zahlenebene mittels komplexer Gréfen fiihrt dabei auf einfache gut iibersehbare Glei- 
chungen. In der I. Mitteilung! ist auf diese Weise das Gelenkviereck behandelt worden. In 
dieser II. Mitteilung soll die algebraische Methode auf die Schubkurbel bzw. Schubschwinge 
angewendet werden. Ferner wird gezeigt, wie diese Methode auch fir schwierigere getriebe- 
synthetische Probleme verwendet werden kann. Als Beispiele sind die Ermittlung der Abmes- 
sungen von Gelenkvierecken aus gegebenen Koppelpunkten und die Getriebesynthese sechsglied- 
riger Kurbelgetriebe bei gegebenen Winkelzuordnungen behandelt. 


2. Zuordnung von Winkeln und Lingen an der Schubkurbel bzw. Schubschwinge. Soll durch 
ein Kurbelgetriebe die Zuordnung von Langen und Winkeln verwirklicht werden®, so kann die 
Schubkurbel bzw. die Schubschwinge verwendet werden. Das algebraische Verfahren mittels 
der Gaufschen Zahlenebene fihrt hierbei zu einer einfachen Darstellung des Problems und 
ermoglicht die Ableitung von Formeln, die einfache Berechnungen erméglichen. Auch Lagen im 
infinitesimalem Abstand, die 
durch Ableitungen dargestellt 
werden kénnen, lassen sich ver- 
wirklichen. 

Wird ein Getriebeglied a 

um die Welle I (Abb. 1) ge- 
dreht, so soll ein anderes Ge- 
triebeglied b gerade verschoben 
werden. Dabeisollennach einer 
gegebenen Funktion den Dreh- 
winkeln @,, g2---p, die Schub- 
langen s,, sg---s, zugeordnet Abb. 1. Schubkurbel. 
sein. Die Koppelgelenke der ; ‘ 
Schubkurbel sind zu ermitteln. Das eine dieser Gelenke ist in seiner Lage festgelegt durch 
und a, das andere durch die senkrechten Koordinaten p und q, die den Vektor t = p +1 q be- 
die Polarkoordinaten stimmen. Werden die Achsen der Gaufschen Zahlenebene durch den Dreh- 
punkt I so hindurchgelegt, daB die reelle Zahlenachse die Richtung der gegebenen Schubbewe- 
gung erhalt, so ergibt sich mit den Bezeichnungen von Abb. 1 fiir die Ausgangslage des Ge- 
triebes folgende Vektorgleichung: 


¢— t+ — 0. (1) 


Dreht sich der Lenker a um den Winkel ¢,, so verschiebt sich der Schieber b um den Weg s,. 
Dann gilt die Gleichung 


gGev—fte*7—t +s, =—0; (20) 
Fiir die zur reellen Zahlenachse spiegelbildliche Anordnung des Getriebes gelten die Gleichungen 
ttre, (2 

und 
Gem enc’ = 5, 0, (4) 


1 Ing.-Arch. 24 (1956) S. 188. 
2 K.-H. Sieker, Konstruktion 6 (1954) S. 351. 


We 
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Die Vereinigung der vier Gleichungen (1) bis (4) fithrt zu der Beziehung 


fi=(—tG—- Nae —« 1256 = tse (5) 
Durch Umformung kann man (5) folgende Form geben 
at(e”—1) +ar(e °—1)as,¢"—as,e +2 p s,s) 30. (6) 


Schreibt man zur Abkiirzung fiir die komplexen Koeffizienten 
e’—1=G(p), & * —-1=3,(9), (7) (8) 
s,0% = B(sp), se” = Gls), (9) (10) 
8; = Gils), (11) 

so erhalt (6) folgende Form: 

at $y) + a-t- $y) —a- FAsy) —a Flsp) +2 ps, — Fels) = 0, (12) 
Sind drei Winkeln ¢, bis y, drei Langen s, bis s,; zugeordnet, so erhalt man aus (12) drei Be- 
stimmungsgleichungen fiir die GréBen a, x, p und q. Eine dieser vier Groen, z. B. « kann frei 


gewahlt werden. Sind vier Winkel-Langen-Zuordnungen gegeben, so erhalt man aus (12) fiir 
y = 1 bis 4 vier Gleichungen zur Bestimmung der vier unbekannten Gréfen. 


a) Zuordnung von drei Winkeln und drei Lingen. Durch eine einfache Umformung 
kann (12) auf folgende Form gebracht werden: 


at lp) +4 FG) + (+t) s, =a Flsy) +4 Flsy) + Bils)- (13) 
Werden fiir g, und s, die zugeordneten Werte ¢, s; Y.S2 und ps, eingesetzt, so erhalt man drei 
homologe Gleichungen, aus denen durch Zusammenfassen folgende Gleichungen gebildet werden 
k6nnen: 


arb = +d, +ad, +s, (14) 
c 1.00 Dye (15) 
(Ce 2) Dia Ged Deeps: (16) 


Hierin sind die 0-GréBen dreistellige Determinanten, die sich aus den gegebenen Gréfen be- 
stimmen lassen. Sie haben folgende Form: 


d=|3,%) BO) s |, (17) 
1 =15(y) Slsp) s, |, (18) 
dD, =1S(y) Blsy) s, |, (19) 
ds =13(y) Sls) s, |, (20) 
6 =| 89) S~) Tse) (21) 
ds =| BY) Bly) Bs) |. (22) 

Durch Zusammenfassen von (14) bis (16) erhalt man nach einigen Umformungen 
a? Gd,—ad?d, +aad,—a2d, + 02d, —aa(d,—d,) ad, +ad3 = 0. (23) 


Setzt man in die Vektorgleichung (23), die ahnlich aufgebaut ist wie (36) in der I. Mitteilung, 
reelle GréBen ein, so ergibt sich entsprechend (46) in der I. Mitteilung 


ce 


a® (d,4 sina + dy4 cosa) —a (2 d,; sina cosa + 2 dy; cos x — dy; + dy. — 2 d,s) 


— (d,3 sina + d,3 cosa) = 0. (24) 

Die Funktion a = f («) nach (24) stellt eine Kreispunktkurve der Bewegung des Lenkers a gegen- 

iiber dem Schieber b dar; sie ist der geometrische Ort aller Drehgelenkpunkte zwischen Lenker a 

sree ae k, die die gegebenen Bedingungen betreffs Zuordnung der Winkel y, und Langen s, 
erfillen. 

Auch fiir den anderen Drehgelenkpunkt der Koppel k zwischen ihr und dem Schieber b kann 

eine Kurve gefunden werden, die sich durch eine Gleichung ahnlich der in (24) darstellen 1aBt. 
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Werden fir einen gewahlten Winkel « die Werte fiir a aus (24) bestimmt, so kénnen die zu- 
gehérigen Werte fiir p und q aus folgender Gleichung ermittelt werden, die sich aus (14) ergibt 
a? b, + a? by + ad, 

ap ; 

Erstes Zahlenbeispiel. Fiir die Zuordnung von Langen und Winkeln an einer Schubkurbel soll 
folgende Funktion zugrunde gelegt werden (Abb. 2): 


t=pt+iq= (25) 


N= C pin ea. (26a) 


Wird der Nullpunkt fiir die Koordinaten y und s an die Stelle c X, =2/6 gelegt und die Kon- 
stante c = | gesetzt, so ergibt sich fiir yy) = 1/4. Mit Beriicksichtigung der Koordinantentrans- 
formation erhalt man dann fiir die Funktion s = f(~) aus (26a), wenn man 


ee aa fe 
setzt,. s -+- - = sin* (y iy 


oder 


und 


CX =Y+CX 


: 1 1 
Ss = sin? (p oc eee 


Werden fiir die Punkte 1, 2, 3 
(Abb. 2) die Winkel 


(26b) 


Ut 


WA MLA 
1 tke P2 6 und = 


eingesetzt, so ergeben sich folgen- 
de Zuordnungen: 


Punkt ¢, ?, Sy 
It 1 
1 15 — — 
ab ae 4 
ney ee as 
6 2 0 IH 
2 
3 45 E + +73) Abb. 2. Darstellung der Funkti ins in? bat ashe 
_ mn b. 2. Darstellung der Funktion y = Csin® ¢ x bzw. s = sin’ (»+2) aaa 


Fir die Berechnung der %-GréBen werden die Sinus- und Cosinus-Werte der Winkel g, und 
ihre Produkte mit den Strecken bend6tigt: 


. 


Punkt sin p, cos gy, s, sin p, 8, COS Y, 
1 0,258 819 0 0,965 925 8 0,064 704 8 0,241 481 5 
2 0,5 0,866 025 4 0,25 0,433 012 7 
3 0,707 106 8 0,707 106 8 0,482 962 9 0,482 962 9 
Aus (7) bis (11) erhalt man dann folgende Zahlenwerte: 
Punkt $,(?) Ty(s9) Sys) 
1 — 0,034 074 2 + 10,258 8190 0,241 481 5 + 10,064 704 8 0,062 5 
2 — 0,133 974 6 + 70,5 0,433 012 7 + 70,25 0,25 
3 — 0,292 893 2 + 10,707 106 8 0,482 962 9 + 10,482 962 9 0,466 506 3 


Die Berechnung der Determinanten (17) bis (22), die mit Hilfe der Beziehungen (64) in der I. Mit- 
teilung durchgefiihrt werden kann, ergibt folgende Werte: 

103d = +7 1,762 078 S, 103d, = + 1,376 73 —i 1,304 133, 

103d, = +. 1,763 974 +i 1,764 261, 103d, = — 2,126 994 +7 2,771 931, 

103d, = — 0,935 754 —i 0,205 263, 103d, = +7 2,935 58. 
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Zur Darstellung der Funktion a = f(a) nach (24) sind folgende Werte berechnet worden: 


nnn Sa LEE EEEESEES EEE! 
2 qd sin x cos A+ 2 an cos? x 
—d d .—0,5d__ 
yt y2 ya 
i 


aa sin Ord, cos p 


ds sin a+d, cos 


ed 


43,44 + 0,549 6 + 0,182 7 0 0 + 0,332 4 
46,68 + 0,354 453 — 0,014 776 + 0,106 822 + 0,528 5 — 0,570 2 
52,50 0 — 0,425 6 + 0,298 3 + 0,700 9 oo 

60 — 0,456 065 — 1,027 520 + 0,540 216 + 0,835 8 ALES 
15 133701 — 2,319 0 + 0,992 3 nO 167 + 0,966 6 
90 — 2,126 994 337-314 + 1,376 73 + 0,733 0 + 0,883 1 


Abb. 3. Darstellung der Kurve a = f(x) nach dem 1. Beispiel. 


Fir den Verlauf der Kurve a = f(x) (Abb. 3) sind folgende GréBen besonders kennzeichnend: Fiir 
& = 52,50° wird d,4 sina + d,4 cos x = 0. Hierfiir ergibt sich der eine Wert fiir a gleich Unend- 
lich, d. h. die Gerade durch den Nullpunkt unter diesem Winkel ist eine Asymptote der Kurve. 
Fir « = 43,44° wird d,3 sina + dy3 cosa = 0 und damit ein Wert von a gleich Null. Da ein 
Wert fiir a nach (23) immer gleich Null ist unabhangig vom Winkel a, ist die Gerade unter diesem 
Winkel eine Tangente der Kurve im Nullpunkt. Der Wert a, = 0,5285 fiir den Winkel « = 46,68° 


entspricht dem im 2. Beispiel errechneten fiir vier gegebene Winkelzuordnungen. 


Als Beispiel ist fiir ~« = 90° und a = 0,8831 aus (25) der Wert fiir t errechnet worden: 
t = — 0,2334 +7 0,694. 
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ee ist p = — 0,2334 und q = + 0,694. In Abb. 4 ist das Getriebe mit diesen Werten darge- 
stellt. 


Um einen Uberblick zu erhalten, wie weitgehend das Getriebe in Abb. 4 die gegebene Funktion 
erfullt, wird der Unterschied 


US Sf (27) 


bestimmt. Hierin ist s, = f(y) die Winkel-Lingen-Funktion des Kurbelgetriebes und s ¢ = f(q) 
die gegebene Funktion (26b). s, kann aus (6) ermittelt werden: 


Sp = A+ FESSOR 28 
Hierin bedeutet 3 y + (28) 
1 ; = : 
eee | Gea) “yp, (29a) 
it Se = 
B= 1 [at (e—1 $ar(e-#— 9), ain 


Abb. 4, Schubschwinge nach dem 1. Beispiel fiir «1 = 90° Abb. 5. Unterschied u der Kurbelgetriebe-Funktion's, = f(y) von 
und a = 0,883. 


der gegebenenFunktion sp=f(9) fiir die Schubschwinge Abb. 4. 


Werden die Gré8en A und B durch reelle Werte ausgedriickt, so erhalt man 


A = a, cos p — ay sing — p, (29b) 
B = (ax p + @y q) (cos p — 1) + (ax q¢ — ay p) sing. (30b) 

Fiir die Schubschwinge in Abb. 4 ergeben sich folgende Zahlenwerte: 
Pe s0e 15) [=27,5)) On ts.) 15 | 225 | 30 |/37,5 | 45-| 52,5 | 60 [ers 
fot = tora 24) 6 | oF iP | Goalies | 0 | +3 | 0 | 20d 255, = 


In Abb. 5 ist die GréBe u in Abhangigkeit vom Winkel ¢ graphisch dargestellt. 
b) Zuordnung von vier Winkeln und vier Langen. Gleichung (12) kann in folgender 
Form geschrieben werden: 


at Sp) + at Fy) —a F.lsy) —AFs (sy) = Gls) —2 psy. (31) 
Setzt man in (31) fiir y, die vier Winkel ¢, bis gy, und fiir s, die zugeordneten Langen s, bis s, ein, 
so erhalt man vier homologe Gleichungen, die in folgender Weise zusammengefaft werden kénnen 


atd=d,+2p ad, (32) 
atdb=d,+2p dz, (33) 
adp=)D,+2p dy, (34) 
ad=bd,+2p dq. (35) 
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Hierin bedeutet 7 
S(sy) Sls)» (36 


( ) 
BAlsy) Fels @)|, (37) 
) 
) 


( 
( 
Sol 


D =|%(9) 
Dy = Oo(s) 


De = | Sy 


Ss 


— = 


Se SS 


ro 


BAlsp) rls) 
Ds = | d(¢) Bs) — Sls) 
Dd, = |S) SY) s TAs q)| - (40) 
Werden (32) durch (34) und (33) durch (35) dividiert, so ergeben sich folgende Gleichungen: 


Ss 


, 


CAI AI ERI_ CAI. CAI 


a d, + 2p dz, Ala 

i Sea ha? (41a) 

= Es (41b) 
ds + 2p dy 


Addiert man die beiden Gleichungen (41a) und (46b) so ergibt sich 


_ ii +2pb by + 2p de (42) 
bap Da bs + 2p by 


Gleichung (42) kann in folgender Weise umgeformt werden: 


(2 p)* Dy Dy + (2 p)? (Ds Dy + D3 Dy — Dy Dy — Dy dy) 
+ (2 P) (D3 else aay cerns yaa 3) — (0, dD, 1D, Dg) == 0. (43) 


Die Koeffizienten dieser kubischen Gleichung sind reell. Die Lésung dieser Gleichung, die z. B. 
mittels des Hornerschen Schemas durchgefiihrt werden kann, fiihrt entweder auf eine oder auf 
drei reelle Wurzeln fiir 2p. Entsprechend gibt es eine oder drei Schubkurbeln, die die vorge- 
gebenen Bedingungen erfillen. Die GréBe q ergibt sich mittels des ermittelten Wertes von p aus 
(41a). Lange a und Ausgangswinkel « der Kurbel erhalt man aus (34). 


Zweites Zahlenbeispiel. Zu den drei Winkeln ¢, bis y3 des ersten Beispiels soll noch ein vierter 
Winkel y, =2/3 hinzugefiigt werden, dem die Lange sy = 3/4 entspricht. Man erhalt also 
folgende . Zuordnungen: 


Punkt °) Py Sy 
ae 1 
1 15 1} 70 
a 1 
aU 
3 45 oy = (12 V3) 
4 60 se Py 
3 4, 


Die Sinus- und Cosinus-Werte der Winkel gy, und ihre Produkte mit den Streckens, haben 
folgende GréBen: 


Punkt sin i 
Py cos Pp, $,° sin D, S$, COS D, 


—eeVO—VOsXKkRFeeele—— 


1 0,258 819 0 0,965 925 8 0,064 704 8 0,241 481 5 
2 0,5 0,866 025 4 0,25 0,433 012 7 
3 0,707 106 8 0,707 106 8 0,482 962 9 0,482 962 9 
4 


0,866 025 4 0,5 0,649 519 1 0,379 
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Fiir die %-Gréfen ergeben sich dann mittels (7) bis (11) folgende Werte: 


Puokt o,(%) &,(s 2) 5, (s) 
1 — 0,034 074 2 + 70,258 8190 0,241 481 5 + 710,064 704 8 0,062 5 
2 — 0,133 9746+ 70,5 0,433 012 7 + 10,25 0,25 
3 — 0,292 893 2 + 10,707 106 8 0,482 962 9 + 10,482 962 9 0,466 506 3 
4 — 0,5 -+ 10,866 025 4 0,375 -+- 10,649 519 1 0,562 5 


Die Berechnung der 0-Determinanten (36) bis (40), die mittels der Beziehungen (88) bis (93) in 
der I. Mitteilung durchgefiihrt werden kann, ergibt folgende Werte: 


10°) = — 489,403, 

10°), = — 150,650 —i 6,4918, 
10°, = + 36,1404 —i 132,602 4, 
10°d; = + 1,621 44—i 335,608 8, 


106d, = — 95,2557 +72 78,5782. 
cot 
‘ool ae 
\ | \ fed BA 7 
\ \ %zo/ \ 
a AY, \ \ 
NIAZ \ BA 
\Wy \ \ \ 
mis a. ~80 
i Ka ats 
e mA \ ‘ 
s roe \ 
S \ \ \ \ 
ie \ . 
eh \ \ 
Ci a hos YZ 
Os 
0,683 
O75 
0939 
~50 
Abb. 6. Schubschwinge nach dem 2. Beispiel. Abb. 7. u = s, = fiir die Schubschwinge Abb. 6. 


Fiir die Koeffizienten der Gleichung (43) ergeben sich folgende Werte: 


1012 D4 Dy = 2. 15.2482, 
1012 (Ds Dy + Dz Dy — Do D4 — Do Du) 5 SIT 
10#2 (Ds D3 —D, Dy — Dy Dg — De D3 — Dy d;) =— 4 166,68, 
1012 (b, Ds + D, Ds) 2 3:868,87 . 


Damit nimmt (43) fir Form an 


(2 p)® — 1,661 02 (2 p)? — 0,273 258 (2 p) — 0,253 727 = 0. 
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Die Auswertung dieser Gleichung, die nach dem H ornerschen Schema vorgenommen werden kann, 
fiihrt zu der Lésung 

2p =+1,87840, p=+ 0,939 20. 
Die beiden anderen Wurzeln sind komplex und deshalb fiir diesen Zweck nicht brauchbar. Aus 


(41) erhalt man 
q = — 0,445 531 


und aus (34) 
a = 0,363 293 +7 0,384 157. 


Hieraus ergibt sich 
a = 0,5281 und « = 46,68°. 


In Abb. 6 ist das Getriebe mit den errechneten Abmessungen dargestellt. 
Der Funktionsunterschied u = s, —s (27) ergibt mittels (26 b) und (28) bis (30) fiir einige 
Winkel folgende Werte: 


0 | 1,5 | 15 oa) 30 | 35. | 45 | 52,5 | 60 | 67,5 


p= | —30 | ae 
105 = tae |==33 | 10) 0 | ea 0 [= 0 ley 0 a 0 eo 


In Abb. 7 ist die GréBe u in Abhangigkeit vom Winkel @ dargestellt. 


3. Anpassung der Ableitungen der Funktionen s = f(@) an die Schubkurbel bzw. Schubschwinge. 
Sollen die Ableitungen der Funktion s, = f(y) an der Schubkurbel den Ableitungen einer gege- 
benen Funktion y = f(x) an einer bestimmten Stelle angepaBt werden, um die Schubkurbel an 
dieser Stelle als Funktionsgetriebe benutzen zu kénnen, so kénnen mittels der Darstellung in der 
Gaufschen Zahlenebene ebenfalls Formeln entwickelt werden, die den fiir die Zuordnung end- 
licher Winkel und Langen entsprechen. Werden die Ableitungen bis zur dritten Ordnung an- 
gepaBt, so ist die freie Wahl einer Gréfe fiir die Bemessung der Schubkurbel méglich. Wird 
dagegen die Ableitung vierter Ordnung auch noch gegeben, so ist eine freie Wahl nicht mehr 
méglich. 

Wird in (12) die Ableitung nach dp bis zur vierten Ordnung gebildet, so ergibt sich, wenn die 
Ableitungen durch Striche gekennzeichnet werden, 


atS(y) +aB'y) —a BF (sp) —ay(sy) +2ps’ —G(s) =0, (44) 
at s(y) +a8'(y) —a F'(sy) —aB'(sy) +2ps” —F(s) =0, (45) 
at e'(p) +a (9) —a W'"(sy) —08'"(sp) +2 ps” — F's) =0, (46) 
at yy) +a R'(p)— 40 B"(sp) —0 B'" (sp) + 2 ps” — F's) =0. (47) 
Da $(p) = e'? — 1 [siehe (7)] ist, erhalt man 
Sy) =ie*, BG) =—e'*, F'(y) =—ie'r, F''"(y) = +e”. (48) bis (51) 
Nach (9) ist 3(sq) = s e'”. Daraus ergibt sich 


a! 


y 
a 


(sp) sells’ +i), oo 
B'(sp) =elv(—s +s" + Bis), Ce) 
B'"(sy) =e (38 48” is 4 B18"), Ce 
B69) =o (6 5" a" his 4 4b"), 55) 


Ferner ist %(s) = s? (11) und damit 
a (s) =2ss8', (56) 


Bs) =2(6 +95", Cy 
a) = 2 (3 s! sl! ae s 3X), (58) 
(8) ex) (3 2 ae 4s’ sf"! ae s suit) f (59) 


Sind fiir die Anpassung drei Ableitungen gegeben, so fiihrt die weitere Entwicklung auf 
Gleichungen, die (23) und (24) entsprechen. Die Determinanten D (17) bis (22) werden aus den 
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Ableitungen der %-GréBen gebildet. So hat z. B. die Determinante D, folgende Form: 


ie”) F@ Fs9) 
m=|8'o) FO Fsy) 


BF") (se) 
Mit dieser Bedeutung der Determinanten hat (23) hier ebenfalls Giltigkeit. 


Ist auch noch die vierte Ableitung s’’”’ gegeben, so entspricht die Entwicklung der fiir die 
Zuordnung von vier Winkeln und vier Langen. Uber die Umformung (31) erhalt man dann zur 
Bestimmung der Grofe p die Gleichung (43), in der die Determinanten ) aus den Ableitungen 
der %-GréBen gebildet werden. So ergibt sich z. B. 


& (s) Sy) Sse) Fs) 
oe (s) Sp) Hse) F(sy) 


d, = z (61) 
os) B() Sse) F's) 


wr") FM FN) Fo) 


Ebenso wie bei dem Gelenkviereck kénnen auch hier Zuordnungen von Winkeln und Langen 
und die Anpassung ihrer Ableitungen kombiniert gegeben und daraus die Abmessungen der Schub- 
kurbel bestimmt werden. 


Drittes Zahlenbeispiel. Die im ersten Beispiel gegebene Funktion 


$= sin® (¢ as 4 == [siehe (26b)] 


soll auch diesem Beispiel zugrunde gelegt werden, um daraus die zu gebenden GréBen zur Ermitt- 
lung der Schubkurbel-Abmessungen zu entnehmen. 

Gegeben sind die Ableitungen im Nullpunkt bis zur dritten Ordnung. Die Ableitungen von 
(26b) fiihrt zu folgenden Beziehungen: 


! . 4 
‘So — 2) sin (r+) cos(p +5] (62) 
s’ =2 eos* (¢ - al — sin? (p + 3}| F (63) 
Dip : a ; We 
s =—8sin(p +4) cos ( + ae (64) 
Fur den Nullpunkt mit g = 0 ergibt sich 
t eine ” we Q 
s'(0) == 3, (0) =1, s”’(0) =— 2/3. 
Werden die Ableitungen in (52) bis (54) und (56) bis (58) eingesetzt, so erhalt man fiir die §-GroKen: 
Abl. &(?) &(s) &(s) 
s As ee 
1s 0O+il oy /3 +70 0 
2. Sy einy 0 eles a 
3. ema inh — {Va +3 | 3/3 


Die Berechnung der )-Determinanten, die in (17) bis (22) dargestellt sind, kann wie im ersten 
Beispiel mit Hilfe von (64) in der I. Mitteilung ausgefiihrt werden. Es ergeben sich folgende 
Werte: 


b = +i37/3 = +1 5,19615, 
 d=+7 510s 
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b= — = }i3 78 


9 5 al 
—_ 


2 4 


— 5,25 +1 5,196 15, 


j3 =—4,5i — 9,093 26, 


Db, = +6+i6/3 =+6 +i 10,392 30, 


db, =—i 6/3 


5= 


22 710:392°305 


Zur Darstellung der Funktion a = f(x) nach Gl. (24) sind folgende Werte berechnet worden: 


2 di sin « cosa+d 1 cos? & 


d_. sina+d_, cosa 
x4 y4 re) ae 
dy tye 0,5d 


d 


eg 8D atd, cos & 


— 8,443 75 
== 455 
+ 0,649 52 


60° 


ys 
60 + 10,392 3 + 18,836 05 
90 +6 + 10,392 3 
120 0 + 1,948 55 
102,03 i. 
120° y 
iv 
Abb. 8. Darstellung der Kurve a = f(«) nach dem 3. Beispiel. 
p =|—30 | —15 |—7,5 


a a, 


a + 0,812 5 
+ 0,866 0 + 0,866 0 
0333 00 


Fir « =116,3 wird 
d,3 sina + d,3 cosa = 0. 

Hierfiir wird ein Wert fir a 
null. Die Gerade durch den 
Nullpunkt unter diesem Win- 
kel tangiert die Kurve im Null- 
punkt. Fiir 0 = 102,38° ergibt 
die Rechnung: a, = 1,1527 
und a, = 0,5842. Der Wert 
fiir ag entspricht dem im vier- 
ten Beispiel errechneten fiir 
die Anpassung der Ableitungen 
bis zur vierten Ordnung. 

Die Kurve a =f (a) ist in 
Abb. 8 dargestellt. Als Bei- 
spiel ist fiir 

“«=—60° und a=-+1] 


aus (25) der Wert fiir ry er- 
rechnet worden: 

t= +1 +7 0,866. 
Damit wird 
Pp=+1 und q= + 0,866: 
In Abb. 9 ist das Getriebe mit 
diesen Abmessungen darge- 
stellt. 

Aus (26b) und (28) bis (30) 
kénnen fur bestimmte Winkel 
gy die Werte fir den Unter- 
schied u = s,—sy (27) be- 
rechnet werden. 


Oilet.5 13 150 \eoolsuliaace | 37,5 | 454 52,5 | 60 | 675 


lou=|+43/+ 2] 0 | o | 0 | +1] +6 [+18] +43] + 86| +157 | +250 | +381 


Diese Werte sind in Abb. 10 in dem mafigebenden Bereich graphisch dargestellt. 


Viertes Zahlenbeispiel. Zu den Ableitungen im Nullpunkt bis zur dritten Ordnung im dritten 
Beispiel soll noch als weitere Bedingung die Ableitung vierter Ordnung zusitzlich gegeben sein. 
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Durch die vierte Ableitung der Strecke s nach d g erhalt man aus (64) 


s’’"' —§ sin (p + a — cos (y -- ol ‘ (65) 
Fir den Nullpunkt mit y = 0 ergibt sich 
Sa So 4 Hl 


Aus diesen vier gegebenen Bedingungen kénnen die Abmessungen des Kurbelgetriebes bestimmt 
werden, wobei es entweder eine oder drei Lésun- 

gen geben kann je nach dem, ob die kubische Wu 

Gl. (43) eine oder drei reelle Wurzeln hat. 


Abb. 9. Schubschwinge nach dem 3. Beispiel fiir = 60° unda = 1. Abb. 10. u = ee: fur die Schubschwinge Abb. 9. 


Aus den gegebenen Werten fiir die Ableitungen erster bis vierter Ordnung im Nullpunkt der 
_gegebenen Funktion erhalt man aus (52) bis (59) fiir die %-GréBen folgende Werte: 


Abl. | 39) 3(s”) 3(s) 
1 0 +71 5 V+ i0 0 
2 —1+4+i0 1+ i/3 = 
3 0—il —7\stis 3/3 
4 130 —10 71013 Sauls 


Die 0-Determinanten (36) bis (40) kénnen wie im zweiten Beispiel mittels der Beziehungen (88) 
bis (93) der I. Mitteilung berechnet werden: 


db =— 432, 
Dee bol | 3 
D, =— 1441 58,5 3, 
bs =— 63 —i153 73, 
pe 99 9 13; 
Fiir die Koeffizienten der Gl. (43) ergeben sich folgende Werte: 
Dy Dy 10044 
De On 00, 0, 0, De == =110:935:, 


DOs 0, Dy) 0, = 0,0, 0, D, = — 10 200, 
Dy.Ds 21D; Ds = +19 845, 
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Gl. (43) erhalt damit folgende Form 
(2 p)’ — 1,088 71 (2 p)? — 1,016 13 (2 p) + 1,975 81 =0. 
Aus dieser Gleichung erhalt man die Lésung 
2p =—1,18277 und damit p = — 0,591 39. 


Die beiden anderen Wurzeln der kubischen Gleichung sind komplex und deshalb fiir diesen Zweck 
unbrauchbar. 
Aus (41) ergibt sich 
q = — 0,528 5 
und aus (34) 
a = — 0,125 22 +7 0,579 74. 
Hieraus folgt 
a = + 0,58433 und o« =102,38°. 


In Abb. 11 ist das Getriebe mit den errechneten Abmessungen dargestellt. Abb. 12 zeigt die 
Funktion u =f() fiir dieses Getriebe nach folgenden aus (26b) und (28) bis (30) errechneten 
Werten: 


Pear 30 | toh arto 


0 | 1,5 | 15 [228 | 30 | 37,5 | 45 | 52,5 | 60 | 67,5 
10°u = | ~104|—38|—9 | 0 | —s |—27|—53 |= 80 | “101 | —113)| tia |= 97 |e 


3. 


0° 
10 


Bild 11. Schubschwinge nach dem 4, Beispiel. Abb. 12. u = Cg fiir die Schubschwinge Abb. 11. 


Fiinftes Zahlenbei spiel. Aus der Funktion 
: Tt 1 
$= sin2 (p tools: TAC 
sind die Ableitungen im Nullpunkt bis zur zweiten Ordnung und damit die Kriimmung der Kurve 
und Punkt 2 (Abb. 2) mit der Zuordnung gy, = 30° und s, = 1/2 gegeben. 
Aus dem ersten und dritten Beispiel ergeben sich folgende Werte fiir die o- GroBen: 


3%) | &(s P) &(s) 

1. Abl. i. 0-Pkt. Oe = V3 +10 0 
2. Abl. i. 0-Pkt. arbre eas = 
Pkt. 2 = V3—14i5 ali ee = 
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Die )-Determinanten der Gl. (17) bis (22) ergeben mittels der (64) I. Mitteilung folgende 
Werte: 


bd = +i (I—5 3) — +7 0,133 975, 

Dd, = + 0,625 — 0,25 73 +4 (2—1,12573) = + 0,191 987 +2 0,051 443, 
D. = — 1,375 + 0,75 ¥3 —i (1 — 0,625 73) =— 0,075 962 +i 0,082 532, 
Ds = —1,125 + 0,625 73 + (0,5 — 0,375 73) = — 0,042 468 —i 0,149 519 , 
Dp 3,5 — 2 30 (2 — 73) = + 0,035 898 +i 0,267 949, 
eee (2065) — — i 0,098 0762. 


Zur Darstellung der Funktion a = f(a) ergeben sich aus (24) folgende Berechnungswerte : 


7 2d__sin&cosa+2d__ cos? 4 F 
a a4 sin aad cos & sie i $s, dis sin a+. cos & a, a, 
DOs + 0,169 128 6 => 0.277 383 3 — 0,113 5877 0,849 3 0,790 8 
60 + 0,165 063 5 + 0,272 116 2 — OFM 538 1 0,885 1 0,763 4 
75,1 + 0,103 614 3 + 0,182 377 8 — 0,079 498 5 0,965 4 0,794 7 
90 + 0,035 898 4 + 0,080 1270 — 0,042 468 3 1,366 0 0,866 0 
92,97 + 0,021 986 7 + 0,060 562 4 — 0,034 675 3 1,942 7 0,811 8 


Die Werte 
Ce a = 0,849, 
== To)", a = 0:795., 
== 97,97, o—= 0,812, 


entsprechen den im sechsten Beispiel errechneten Werten fiir vier gegebene Bedingungen. 

Fir den Winkel « = 97,6° ist d,4 sina + dy4 cosa = 0. In dieser Richtung lauft die Kurve 
ins Unendliche; die Gerade durch den Nullpunkt in dieser Richtung ist eine Asymptote der 
‘Kurve. 

Fir den Winkel « = 99° ist d,3 sina + dy3 cos « = 0; die Gerade durch den Nullpunkt in dieser 
Richtung ist Tangente der Kurve im Nullpunkt. 

In Abb. 13 ist die Funktion a = f(x) dargestellt. Als Beispiel ist fiir « = 90° und a = 1,366 
aus (25) der Wert fiir r errechnet worden: 


| feo 7 1183.25 


Damit ist p = 0 und q = + 1,183 2. Abb. 14 zeigt das Getriebe mit diesen Abmessungen. In 
Abb. 15 ist die Funktion u = f(g) fiir dieses Getriebe dargestellt nach folgenden mit (26b) und 
(28) bis (30) berechneten Werten: 


30 8%)5 | 45 1152.5 | 60 | 67,5 


p= | —30 ea j= 


102 u = | 250) = ie 0 


Ope ageen is | 22,5. | 

et te 
_ Sechstes Zahlenbeispiel. Zu den Bedingungen des fiinften Beispiels ist als weitere noch die erste 
Ableitung im Punkt 2 (Abb. 2) und damit die Steigung der Kurve an dieser Stelle gegeben. 


Fir Punkt 2 mit g, = 7/6 ergibt sich aus (62): 
/ 1 a 
S, = Gp y3 . 
‘Die Kurve hat also im Punkt 2 dieselbe Steigung wie im Nullpunkt. 


Die Getriebeabmessungen kénnen wie im zweiten und vierten Beispiel aus (43), (41 a) und (34) 
berechnet werden. 


| 
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Abb. 14. Schubschwinge nach dem 5. Beispiel fiir ~« = 90° 
und a = 1,366. 


58,9° 


Abb. 13. Darstellung der Kurve a = f(x) 
nach dem 5. Beispiel. 


Abb. 15. u = 3 f fiir die Schubschwinge Abb. 14. 
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Die Werte fiir die @y-GréBen kénnen dem zweiten und vierten Beispiel entnommen werden: 


—X—“"_—Sooeeeeeeeeeeee—————— — — 


3(%) &(s ) 

1. Abl. i. 0-Pkt. Coal + V3 + i0 

2. Abl. i. 0-Pkt. ie 0 147)3 

i. al 1 y= 1 

Pkt. 2 —/3 — —- — | — 

5 V3 Meas Tira, 
: es be 1 li 

1, Abl. i. Pkt. 2 Se pe re 
1 3 +i 3 \3 2 +1 5 3 


&(s) 


Ale bo| oo 


1 y= 
V3 


Die Berechnung der )-Determinanten nach (36) bis (40) fiihren mit Hilfe von (88) bis (93) der 


I. Mitteilung auf folgende Werte: 


b => (—21 +1273), 103) = — 53,8475, 
db, = (8 — 273) +12 (26 —15 73), 10ehi ee 

b, = | (—7 + 473) —i + (33 — 1973), 105) 

ds = 5 (45 — 26 73) i (9 — 53), 103 0 

D, == (—59 + 34/3) +4 (26— 1573), 10°, = 


Die Koeffizienten der Gl. (43) ergeben sich wie folgt: 
d, dy 


Ds Dy + Ds Dy — dy Dy — Dy Dg 


29,0064 +7 3,60713, 
— 8,974 60 —z 11,379 4, 
— 4,165 13 —i 42,468 3, 


— 13,7841 ++i 2,404.75. 


1 ax 
= +, (8 300 — 4792/2), 


1 — 
— 55 (4 874 —2 8143), 


i y=; Dy — 0, B;— Da Dy — 0a dg =F -, (627 — 3023), 


D, Dp +d, ds = +2 (362 — 209 y3). 


Gl. (43) erhalt mit diesen Zahlenwerten folgende Form: 
(2 p)? — 1,437 (2 p)? — 0,189 7 (2 p) + 0,325 = 0. 
Diese kubische Gleichung hat drei reelle Wurzeln 
2p, = +0,495, 2p,—+1,408, 2p, =— 0,466. 
Damit wird 
Pi = + 0,248, Po = + 0,704, Ps = — 0,233. 
Aus (41a) und (34) ergeben sich die entsprechenden Werte q, a und a zu 
g1 = + 0,7445,. gg = +0,5408, 9g, = + 90,5574, 
a, = + 0,7933, a,=—+0,8476, a, = + 0,8106, 


a =75,1°, Ot, = 58,9°, OOP 


In Abb. 16, 17, und 19 sind die drei Getriebe dargestellt. Das Getriebe in Abb. 16 ist praktisch 
nicht brauchbar, weil sich die Getriebeglieder nicht von der einen Stellung in die andere ohne Aus- 


einanderbau iiberfiihren lassen. 
Fiir die beiden anderen Lésungen ist die Funktion u = f(g) in Abb. 18 


und 20 dargestellt. 


Den Kurven liegen folgende mit (26b) und (28) bis (30) berechneten Werte zugrunde. 
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p= |—30|—15|—75| 0 eaiee | 22,5 | 30 37,5 | 45. | 52,5 | 60 | 67,5 

Ce ere eth 4) pet rom 6 +2 |+13| +44 | + 104] + 198 
Lésung 3: 

g= |—30|—15|—%5| 0 |, 4,5 || 15 | 22,5 | 30 37,5 | 45 | 52,5 | 60 | 67,5 

10° u = | ely] eae == | 0 | 0 0 | 0 | 0 0 |+2|+43|—7 ante 


om 0,704 = 


Abb. 16. Schubschwinge nach dem 6, Beispiel Lésung 1. Abb. 17. Schubschwinge nach dem 6. Beispiel Lésung 2. 


10° 


495, 60 


Abb. 18. u = So ‘f fiir die Schubschwinge Abb. 17. Abb. 19. Schubschwinge nach dem 6. Beispiel Lésung 3. 


4. Ermitthung der Abmessungen von Gelenkvierecken aus gegebenen Punkten der Koppel- 
kurve und Anwendung auf die Synthese sechsgliedriger Kurbelgetriebe. Wie aus analytischen 
Untersuchungen bekannt ist, ist die Koppelkurve eines viergliedrigen Kurbelgetriebes eine Kurve 
sechsten Grades!. Eine solche Kurve ist durch neun Punkte in einer Ebene bestimmt d. h.: 
Werden neun Punkte in einer Ebene gegeben, so muB sich durch diese eine ganz bestimmte 
Koppelkurve hindurchlegen lassen. Ferner miissen sich aus diesen neun Punkten drei Gelenk- 


1 H, Alt, Z. angew. Math. Mech. 3 (1923) S. 13. 


he 2 
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vierecke nach Lage und Abmessungen ermitteln lassen, deren einer Koppelpunkt die Koppel- 
kurve beschreibt; denn nach Roberts! la&t sich jede Koppelkurve durch drei Gelenkvierecke ver- 
wirklichen. Das an sich bekannte Problem ist bisher weder algebraisch noch geometrisch gelést 
worden. Es ]a8t sich jedoch mit der komplexen Darstellung in der Gaufschen Zahlenebene lésen 
und fiihrt zu verhaltnismafig einfachen Beziehungen. 


Das Problem soll folgendermafen formuliert werden: Die Punkte C, bis C, (Abb. 21) sind in 
der Gaufschen Zahlenebene gegeben durch die Vektoren tr, bis ry. Im allgemeinen Ansatz soll y als 
Index verwendet werden. Zur Bestimmung des Gelenkviereckes A, ABB, mit dem Koppelpunkt 
C, dessen Koppelkurve durch die Punkte C, bis C, hindurchgehen soll, sind neun Bestimmungs- 
stiicke erforderlich: Die Drehpunkte im Gestell 4, und By sind durch die beiden Vektoreny, und 
1, erfabt, die je zwei Bestimmungstiicke, nimlich GréBe und Richtung, enthalten. Die weiteren 
Bestimmungsstiicke des Gelenkvereckes sind die Lenkerlangen a und b, die Langen c und d des 
Koppeldreieckes 4 BC und der Winkel (6 — y) zwischen den Langen c und d. 


—50 


Abb. 20. u = ies *f fiir die Schubschwinge Abb. 19. Abb. 21. Gelenkviereck mit Koppelkurve des Koppelpunktes C. 


Mit den in Abb. 21 eingetragenen Bezeichnungen kénnen folgende Vektorgleichungen auf- 
gestellt werden: 


a ev =T, —'%—C€ ely t%) ? (66) 
be’ — t, —t, —d eilotxy) ; (67) 


Denkt man sich das Getriebe an der reellen Zahlenachse gespiegelt, so erhalt man entsprechende 
Gleichungen mit den konjugiert-komplexen Gréfen: 


em the Cem). (68) 
a a ae (69) 


Werden jeweils die Gleichungen (66) und (67) und die Gleichungen (68) und (69) miteinander multi- 
pliziert, so ergibt sich i s is 
a? = (iy a Chan ei ltt *y)) (ty a gt aan Cre. ”) ? (70) 


ee eee ON) (Cp de ee (71) 


Um aus diesen beiden Gleichungen den Winkel x, eliminieren zu kénnen, werden sie durch Um- 
formung und Zusammenfassung auf folgende Form gebracht: 


et) (9 BAB) =—q- A—pB, (72) 
e+) (Of B— AB) =p B— GU. (73) 


1 §. Roberts, Math. Soc. Proc. 7 (1875) S. 14. rs 
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Multipliziert man diese beiden Gleichungen miteinander, so erhalt man 
(AB — WB)? + (¢AX—p B) (qU—p B) =0. (74) 


Hierin bedeutet, wenn man fiir c e'(—”) = ¢ setat; 


=i) ce ee (ee (75) 
S=(t,—t:)d, 6 =(t,—t,)d, (76) 
p = (t, —%a) (ty —Ta) + 2? —a?, (77) 
Gq (Ep te) (te (78) 


Setzt man in (73) y = 1 bis 9 ein, so erhalt man neun homologe Gleichungen zur Bestimmung 


der neun Unbekannten rq, Q, rp, y, a, 6, c, d, d — y. 

Werden die Klammerausdriicke in (73) unter Beriicksichtigung der Beziehungen (75) bis (78) 
ausmultipliziert — eine etwas miihsame, wenn auch nicht schwierige Arbeit — und die Glieder 
geordnet, so erhalt man eine Gleichung von folgender Form: 


Ey) (6,71, +6, 1, + Ey +O. +@, 733 +, ft, 4 6.7%) Or eae 
— Fit — at — he ts — Oat, — a t, — FO (79) 


Die &- und %-GréBen sind Funktionen der neun Unbekannten, von denen fiir die weitere Ent- 
wicklung besonders die €-Gréfen interessieren. Zur Abkiirzung ist in den folgenden Beziehungen 


rteft&—@=u?, (80) 
7? + d? —b? = v? (81) 

gesetzt: 
E, =(c—d)(c—d), (82) 
€, =—(c—d)[(c—4) (1. + ts) + (C—T. a), (83) 


E, = @ (c—t)? + (€ —d) (co? —‘du®) + (¢ —d) (¢ »? — du?) + (t, € — tod) (Tt. C—% 4), 
+ (¢ — 4d) [(ta + ¥s) (te C— Ta d) + (Ta Ts C— Ta Vs A)] 


+ (c —D) [(f. + Ts) (t2¢ —t.d) + (ta tC — 1. td), (84) | 
€, =(C —Dd) [(%. + ts) (tC — Tad) +1. (¢ —d)], (85) 
GC, = — 2 d#(¢ e)(t; c—1,, c) —(€ v2 — du?) ce —_d) (S, + ty) 4 tc — tea) 

— T, t,(¢c — d) t4 ¢ — tq d) — (Ty C— Td) (Ta Uy C — Va Tp A) 

— (Ta + Ty) (C — d) (ta tT, C—T, yd), (86) 
C, = —T, t (¢ — @) (% ¢— Td). (87) 


Werden in (79) fiir 1, die gegebenen Vektoren 1, bis ty eingesetzt, so erhalt man neun homologe 
Gleichungen. Die %-GréBen, deren Zusammensetzung aus den neun Unbekannten hier nicht 
interessiert und deshalb nicht angegeben zu werden braucht, lassen sich durch Zusammenfassen 
von jeweils sechs der neun Gleichungen eliminieren. Zu diesem Zweck werden die %-Gré8en 
auf die linke und die €-Gré®en auf die rechte Seite der Gleichung gebracht. Gleichung (79) 
erhalt damit folgende Form: 


Fir} + Jv; ae + st, + Ft, + F, = Er + Grit, +, rT, 


+£yr7+G&hu+Grt+G, rt, +671, +1 + G3. (88) 
Wird y = 1 bis 6 eingesetzt, so erhalt man durch Zusammenfassen der sechs Gleichungen 
F, bd, = £, d¢, + ©, b%, + Eds, + Ey ds, +... & d¢,, (89) 
20 = ED ee eee (90) 
Be db) = E, di, + Gd +, (91) 
O01 = Ed! + G ds +... , (92) 
Bs di = E, di, + Gd, +e, (93) 
F,d,=£,df,+ © df +..-. (94) 
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Werden sechs von den neun Gleichungen in anderer Weise zusammengefabt, z. B. fiir y = 1, 2, 3, 
4,5, 7, so erhalt man 


F,d, = E,d¢, + ©, ds, + Gy d¢, + E,d2, + --- dz, (95) 
%e de = E,d', + E, d?,, | (96) 
fis De == Bde ess, (97) 
Peta = Os Oe ss, (98) 
Ogee HDS doses, (99) 
Fe 0,=— EH, if, ae +: . (100) 


Mit anderen Kombinationen der y-Ziffern lassen sich andere Gruppen von sechs Gleichungen 


bilden. 

Wird die erste Gleichung der ersten Gruppe durch die erste Gleichung der zweiten Gruppe 
dividiert, so fallt F', heraus und es ergibt sich nach einigen Umformungen 
E, (0¢, Ds — D4, D1) + E, (04, D. — D4, D1) + Ee (DF, Dg — dS, D1) +++ Ge (D3, d2—bs,d,)=0. (101) 


Werden auch die zweite, dritte bis sechste Gleichung der beiden Gruppen in derselben Weise 
zusammengefaBt und umgeformt, so erhalt man 


E, (d?, De r- OF; D,) si C, (di, De ae On D;) hase Gy (Oy, De aa, Des i) =0 (102) 
E, os De —y Das D4) Seon (103) 
usw. 


Nimmt man eine dritte Gruppe dazu, so kénnen weitere Gleichungen von folgender Form gebildet 
werden: 


E, (df, Ds i Ds, oy) sie C, (df, Ds a Dee D1) Ae a G, (df, Ds pas dso d,) (104) 
E, (b°, D3 — 08, D,) +--- (105) 
usw. 
Auf diese Weise kénnen also leicht neun Gleichungen gebildet werden von folgender Form: 
E, Ay + © dy + & Aig + +++ © Ai =9, (106) 
Ey Any + Gy Ase + Gy Aog + +++ Gs doy = 9, (107) 
: (108) bis (113) 
Ey Agi + ©z Age + Ee Ag, + +++ €s oo = 9, (114) 


Werden die Glieder einer Spalte, z. B. der ersten, auf die rechte Seite gebracht, und die neun 
Gleichungen zusammgengefaBt, so erhalt man 


6D =—F,Dodr6,-—-hY G=—-H, (115) 
E,=—E, = (116) 
= €&=—-H%, (117) 
6=--Eh2 €=—£ =, (118) 
6=-—k, =  G=—-h2. (119) 


Wie man zeigen kann, sind die GréBen D,/D, D,/D, D,/D und D,/D jeweils konjugiert-kom- 
plex zu D,/D, D4/D, De/D und D;/D. Es werden deshalb folgende abgekiirzte Bezeichnungen 


eingefihrt: 


t 


2 a a Dt gg, 22 = Gy usw. 120 
ie Og EY (120) 


18* 
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Die Gleichungen (115) bis (119) nehmen damit die Form an 


¢, = — EF, 41; €, =— E, 4,, (121) 
E, =— E, q2, (122) 
G, =— E, 43, ©, =— E,4;. (123) 
€,=—E,q, €, =— E,4,, (124) 
€, =— £, 45; €, =— E, qs. (125) 


Werden die Werte fiir ,, &, und G, nach (83), (85) und (87) in (121), (123) und (125) eingesetzt, 


so erhalt man 


(c—d@[C—d) (te +0) + (ut —td)) =F, (c—ade—4), (126) 
—(¢—) [( —d) ta ty + (te + ts) (tC — Ta d) = Ga (C4) (CA), (127) 
(C2) tan (ae ==) = tp (C= aCe (128) 


Wird aus diesen drei Gleichungen durch Zusammenfassen der Ausdruck t, ¢ —1t, d eliminiert 
und durch Abkirzug t, + t, = © und t, t, = § gesetzt, so erhalt man 


ieee fag a eee (129) 
c—d S ss 


Durch Umformung ergibt sich 


B==* = 6 (G — 4) —43, (130) 


— 6 
COP GS (0) — GSE (Gris dg) 0 (131) 


Die kubische Gleichung ergibt drei Lésungen und damit drei Gelenkvierecke, die die gestellten 
Bedingungen erfiillen, wie es nach Roberts auch sein muf. 


Ist aus (131) und (130) die Summe und das Produkt der beiden Vektorenr, und 1, berechnet, 


so kann auch jeder Vektor einzeln bestimmt werden. 


Aus (129) ergibt sich die Beziehung 
_ ta t= 4s (132) 


Sind aus (130), (131) und (132) die GréBen r,, tv, und ¢/d berechnet, so sind auch die konjugiert- 
komplexen GréBen f,, t, und c/d bestimmt. Sie lassen sich auch aus (121), (123) und (125) nach 
demselben Entwicklungsgang ermitteln. 


Zur Bestimmung der Groen u, v und c bzw. ¢ werden die tibrigen noch nicht benutzen von 
den neun Gleichungen (122) und (124) herangezogen. Sie kénnen wie folgt geschrieben werden 


#02) + G9 Ge—#) + (pa) Gow) + (ual —8) 
+(§-a)[@. +0 (25-2) +(.52—r4)] 


ae Ge 1 Ce ts) a) a ( te — Te i)| => tee 1) eo 1) » (133) 


set —2ee 5) + Gee) E—a)ee+9 + (65-2) 


+ (Te +%)($—1) ae i) = a(S =a (5 —1}; (3% 
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se(E =i) 2s (ol fg— apne) 
+ Veale Eostig es +(6$—1) (tt —7, i] 


t (ta 4 wo) ($—1)(u$—1 8) =m (£—1)(E— 1). (135) 


Aus (133), (134) und (135) kénnen die Ausdriicke c2, “v2 
mit denen dann auch u, v und ¢ bekannt sind. g 
Von der Durchrechnung eines Zahlenbeispieles mute Abstand genommen werden. Wenn 
man sich den Gang der Rechnung vergegenwartig, so erkennt man, daB die Rechnung zwar 
nicht schwierig ist, weil nach iibersichtlichen Schemen gerechnet werden kann, aber doch erheb- 
lichen Zeitaufwand fiir den einzelnen Bearbeiter bedeuten wiirde, wenn keine besonderen Hilfs- 
mittel zur Verfiigung stehen. Die aus den gegebenen t,-GréBen zu berechnenden q-GréBen sind 
Quotienten neunstelliger Determinanten, deren Glieder sich aus sechsstelligen Determinanten 
zusammensetzen. Die Rechnung kann also nur in einem ertraglichen Zeitaufwand durchgefwhrt 
werden, wenn ein Stab von Bearbeitern zur Verfiigung steht oder mechanische oder elektrische 
Hilfsmittel etwa in Form eines ,,Determinantenrechners“* verwendet werden kénnen. Dieser 
wie auch der nachste Absatz ist deshalb nur als Anfang dafiir gedacht, in weicher Weise sich die 
algebraische Methode auch auf kompliziertere 
getriebliche Anordnungen mit Vorteil anwenden 
lassen, ohne den Anspruch auf erschépfende Be- 
handlung des Problems zu erheben. Er laBt ein 
weites Feld fiir weitere Untersuchungen und 
Forschungen offen. 


¢ ° 
— u* und ri v* —u* bestimmt werden, 


Ze 


Abb. 22. Sechsgliedriges Kurbelgetriebe in der Anordnung Abb. 23. Kurvengetriebe als Ersatzgetriebe fiir das sechsgliedrige 
eines Zweistandgetriebes. Kurbelgetriebe in Abb. 22. 


Hier soll nur auf eine getriebesynthetische Anwendung des in diesem Abschnitt behandelten 
Problems hingewiesen werden: Bei der Bewegung eines sechsgliedrigen Kurbelgetriebes in der 
Aufstellung eines Zweistandgetriebes (Abb. 22) bescbreibt der Kurbelzapfen A des Gliedes a 
relativ zu dem anderen im Gestell drehbar gefiihrten Glied b eine Kurve. Wird das Getriebe in 
kinematischer Umkehr auf das Glied b als Gestell aufgestellt gedacht, so ist die Kurve die Koppel- 
kurve eines Gelenkviereckes, das auBer b aus den Gliedern g, f und c gebildet wird. Beim Antrieb 
von Glied a wird Glied b im sechsgliedrigen Kurbelgetriebe ebenso bewegt wie in einem Kurven- 
getriebe nach Abb. 23, einem Ersatzgetriebe fiir das Kurbelgetriebe *, wenn die Kurve nach 
der Koppelkurve von A relativ zu b gestaltet wird. 

Sollen nun acht Winkelstellungen y, des Gliedes a acht Winkelstellungen y, des Gliedes b 
zugeordnet werden, so lassen sich bei freier Wahl der Gliedlange a und seiner Anfangsstellung « 
aus den gegebenen Winkeln gy, und y, neun Punkte der Koppelkurve festlegen, die Punkt A 
relativ zu b beschreibt (Abb. 22). Mit Hilfe der in diesem Abschnitt entwickelten Beziehungen 
kénnen dann die iibrigen Abmessungen des sechsgliedrigen Kurbelgetriebes, die sich auf das 
Gelenkviereck b, f, c, g beziehen, ermittelt werden. Dabei werden die Langen auf die Gestell- 
lange als Einheit bezogen. Eine Funktion py = f(y) kann also mit einem sechsgliedrigen Kurbel- 
getriebe in neun Punkten verwirklicht werden. 

5. Winkelzuordnungen am sechsgliedrigen Kurbelgetriebe. Mit der Zuordnung von acht 
Winkelpaaren ist am sechgliedrigen Kurbelgetriebe die héchstmégliche Zahl der Zuordnungen 


1 K. Hain, Feinmech. Praz. 49 (1941) S. 219. 
2 G. Kiper, Konstruktion 6 (1954) S. 300. 
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nicht erreicht, weil zwei GréBen a und « noch frei gewablt worden sind. Wird diese Freiheit dazu 
ausgenutzt, weitere Winkellagen zu verwirklichen, so kénnen zehn Winkelzuordnungen (vorge- 
schrieben werden, aus denen dann die Abmessungen des Kurbelgetriebes, die sich durch zehn 
BestimmungsgréBen festlegen lassen, dargestellt werden kénnen. Damit kénnen von einer 
Funktion y =f(q) elf Punkte verwirklicht werden. Fiir die zehn das Kurbelgetriebe festlegenden 
GréBen werden gewahlt (Abb. 22): a, a, by, Bi, be, Bos C1» V'19 C29 Yor Diese zehn GréBen lassen sich 
als Vektoren durch folgende komplexen Groen darstellen: a, 61, bg, C1, C2. Die Lésung dieses 
Problems wird hier am Zweistandgetriebe gezeigt, kann aber sinngemaB auch auf die anderen 
Formen der sechsgliedrigen Kurbelgetriebe angewendet werden. 

Wie an Abb. 22 verfolgt werden kann, kénnen mit den im Bild eingetragenen Bezeichnungen 
folgende beiden Vektorgleichungen aufgestellt werden: 


ae’ rand e Hefei era =0, (136) 
eA Ce ee be te (137) 


Die Gleichungen gelten fiir eine Stellung des Getriebes, in der sich Glied a um den Winkel 9, 
gegeniiber der in Abb. 22 gezeichneten Stellung weitergedreht hat. Das Glied b hat sich dann 
um den Winkel y, und das Glied c um den Winkel &, gegeniiber der gezeichneten Lage gedreht. 
Die Winkel ¢, und 7, fiir die Koppeln f und g werden gegen die Waagerechte gemessen. 

Denkt man sich das Getriebe an der reellen Zahlenachse gespiegelt, so gelten entsprechende 
Gleichungen fiir die konjugiert-komplexen GréBen 


Geutcyes uetfe (bnew lee (138) 
ae” + Cy eit =P ei =f, ges (emesis (139) 


Werden die Gleichungen (136) bis (139) nach den f- und g-Gliedern aufgelést und dann (136) mit 
(138) und (137) mit (139) multipliziert, so ergibt sich 


fila (lee bees e ee —tnee (Lee bieo en — 10, cae (140) 
eg? = (14+ 3b, e = qe — Ce e-”) (oe bs ea Une Ven a Ce e*) ; (141) 


Diese beiden Gleichungen kénnen durch Umformung zwecks Eliminierung der GréfRe é, auf 
folgende Form gebracht werden 


a eae Rea (142) 
Ge Ue +t,e " U,—=B,. (143) 

Werden diese beiden Gleichungen so zusammengefaBt, daB &, herausfallt, so erhalt man 
(C1 Co Wi We — Cy Cy Wy ,) + (¢, U; B, — ¢; Us B,) (¢, % B,—c,%, Bi) =0. (144) 


Hierin bedeutet 


0S 1 hye ee, (145) l= 1 be %—ae (146) 
US 1e hres 2 —aens (147) Yj =1+hes"—aen”, (148) 
Hepes 7 SES eS peo (149) B, = ty oC (150) 
wy = bt +o +e —f2—1, (151) u, = 62? +a 4+ c2—g?—1, (152) 
C, =a 5, ea Deg untae) CG, =a b,c "eG bee 9 ad) 


Setzt man fiir die Anfangslage @,, y, und &, gleich Null, so erhalt man aus (140) und (141) fiir f 
und g die weiteren Beziehungen 


PS (1b 06) (eo, a ee (155) 
g? =(1+6,—a—c,)(1+6,—a—é). (156) 
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Werden die Klammern in (144) durch Ausmultiplizieren unter Beriicksichtiung von (145) bis 
(154) aufgelést, so erhalt man nach dem Ordnen der Glieder eine Gleichung folgender Form: 


le 


+ ©, 3 Fy—¥y) 


i G, Ace EAC oem) mm ©, giety) G, ga i(28, + ¥y) 
6, ei Prt 2) L &, ey +2 ey) ©, 0% ae G, on 8% 
. err) G, es My) obi 6, eit? ae 8% 
—- %e ee) es % eM 2y,) Xs ei2%—3 vy) Ss gia —3 9) 
= BS See UG ade Xe ety — Fy) % eil(%y—3y%) __ 5 oi 3%) (157) 
best Gy en Wy) <<: G Fah ame) £5 G, etl — 2%) es 6, il %—2 vy) 
a ®, el? ty) pan 6. ei Fwy) Cn ei Py — Vy) ny G, eh ,) ue 
ot Ge?” a D eit =O. ei2My ee Se eit My 
eae Gee er 9. en ey a8 ge? pad Sa oni Eine 


Die C-, %-, G- und -GréKen sind Funktionen der unbekannten fiinf Vektoren a, DysiUns Cisa las 
Sie enthalten auBerdem die Koppellangen f und g, in den Gré®en u, und uy, die aber nach (151) 
und (152) aus den fiinf Vektoren darstellbar sind. Fir die &-GréRen, die fiir die weitere Ent- 
wicklung besonders intressieren, gelten folgende Beziehungen: 


¢, =— a3 by be (Cy — Ce) (Cy bs — Ce b,) > (158) 
©, = + a? (Cy — Ce) [( C1 — Cz) (1 + by) + (C1 Be — Ce Bi), (159) 
©; = — a (Cy — Ce) [(C, — Ce) By be + (6, + Da) (C1 be — Ce bi], (160) 
©, = + by by (C, — Ce) (C1 be — 2 by), (161) 


GC; = + a? [(Cy Co by — Cy Cy Bg)? + 2 (Cy by — Ce be) (€; by — Ey by) +4, By (Cy —C2) (Co — Gi) 
+ (1 + a*) (¢y by — Co By) (€, by —Ty By) + 4, By (¢: — 0) (€ By be — Ca By ba) 
+ 2 (6; +b.) (cy — C2) (C, by —Ce by) + bi b, (C, — Cz) (Cy Up — Ce Uy) 
4 (Ce by — Ce 6) (¢, by wu, — ¢, by u;) 
C, = — a? (Cy — Ce) (C, — Co)- 


(162) 
(163) 


Setzt man in (157) y = 1 bis 10 ein, so erhalt man zehn homologe Gleichungen. Durch wieder- 
holte Anwendung der im vorigen Abschnitt benutzten Methode [siehe (88) bis (105)] kénnen die 
§-, G- und %-GréBen eliminiert werden, so daB sich zehn Gleichungen von folgender Form 


ergeben 


€,4,+6, Ai2+ 6, Ais + &, Ay, + ©s Ais +6; Aig + &4 Ay, 

+6, Aig + C5 Ay +65 Aro =—&, A), 
€, 4+, Asa + ©, Ao3 +E, Asst °*: 
€,Ant+G, dat 


(164) bis (173) 


Werden diese zehn Gleichungen zusammengefaBt, so erhalt man 


6, A+ Aare: 
~ D % 
¢, = — Cs =— Ge» (174) 
&, —a &, Qe > (176) 
On i &, Qs > (178) 
Cy = — Ce Qa» (180) 
¢, = oo &, q5 2 (182) 


€, = 
¢; = 


— 6, qi > (175) 


) 
ed OF qs f (179) 
— ,q,,- (181) 
— €,q,- (183) 
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Setzt man die Werte von (158) bis (161) und (163) in (175), (176), (178) und (180) ein, so ergibt 


sich 


— a b, be (Cy — Ce) (C Be — Ce by) = a? (¢, — Cp) (Cy — C) > Gy» , (184) 
+ a? (€; — ©) [(¢, — Ce) (By + Be) + (Cr ba — C2 by)] = a5 (Cy — C2) (Cy — Ce) qa» (185) 
— a (© — Ce) [(¢r — Ce) by be + (Or + Ba) (Cy Bz — Cy by)] = 8 (C1 — Cy) (Cy — Ce) Ga» (186) | 
+ b, by (Cy — Ce) (1 By — Ce By) = 03 (¢, — Cy) (Cy — Ce) Ga. (187) 
Die Zusammenfassung von (184) und (187) fihrt zu der Beziehung 
a? = — a (188) 


* Werden (185), (186) und (187) zusammengefaBt, so ergibt sich 


Cy bg— Cy by __ 5 — _ gst hy be a as 
so aerated ea L (by 4 be) = Ete bee bh (189) 


Setzt man zur Abkiirzung b, + 6, = © und b, b,= §, so lassen sich die Beziehungen in folgender 


Form schreiben 


neo. 
eg ce (190) 
SG? — 26? a q, + Ga? (q5 — ga) + 0° (Ge G3 + qa) = 9, (191 
¢ bib, — oa, 
Cy =a Ja (192) 


Gleichung (182) kann unter Beriicksichtigung von (162) wie folgt umgeformt werden 


ot 225) +2(6,— : ba) (= be i) A5, bs ( —2\(2—1) 
1 . ; 


52,155) te alle bo 6} =| Dabs (tere ~—1) “ (2 b.—,) (by ma — 7 bau) 


Eine entsprechende Umformung kann mit (183) unter Beriicksichtigung von (162) vorgenom- 
men werden: 


c, (2 2 by — be) 42 @ by —b,) (be — 2 by) Soh b(1—2)(2 — 7 


C3 CANN 


@ 2 Gi 2 2 ¢, | 
+2 (6; + 6)(2— )(6. ; by) 4 6, b.(2 1) (4 — 2m) + (be —26))(2 b, ty bi) | 


pa ok _ Ca) é , | 
ae e alle 1 ee (194) 


Aus (151) und (152) erhalt man unter Beriicksichtigung von (155) und (156) fiir u, und u, die 


Ausdriicke 


u,=¢(1+5, a) + Cy 2 (1 + b, —a) —(1 +6, —a)—(1 + 6, ~~) +. 0b, +45, (195) 


¢ 


rate ars een =) a 
= Tee t+ by — a) r C2 (1 + b,—a) —(1 +b,—a)—(1 +b, —@) +06, +45,. (196) 
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Werden diese Ausdriicke in (193) und (194) eingesetzt, so kénnen diese auf folgende Form ge- 
bracht werden: 


Gilg ts Caly —= Crs Ue 1D, (197) 
Gi Veuste Cato =e Gy) Cy Does 5 (198) 
Durch Zusammenfassen dieser beiden Gleichungen ergibt sich 
Cy (Jy tp —Jo ty) = Cy Ce (dy tp — dy t,) + (tv,t, — wy t,). (199) 
Co (la te — Jo ty) = Cy Ce (fy D2 —Jo dy) + (fy 1, — J 10;)- (200) 


Durch Multiplikation der beiden Gleichungen (198) und (199) miteinander erhalt man 

(C1 C2)? Dy tg — Vp th) (fr De —Je 0) 

+ Cy Cy [(Dy ty — By ty) (fy a —fo 1) + (1, ty — Wy ty) (f4 Yo —Jo 1.) — (fs te —Jo t,)7I 

+ (tv, tf, — te t,) (J, 2 — fy.) = 0. (201) 


Aus dieser Gleichung kann ¢, ¢, und mit Hilfe von (199) und (200) dann c, und ¢, ermittelt werden. 


Der Rechnungsgang ist demnach folgendermaBen: Die GriBen gq, bis q, und q, bis q, werden 
aus den gegebenen Winkelzuordnungen berechnet. Aus (188) erhalt man @ und damit auch a 
und zwar im allgemeinen drei Werte. Aus der kubischen Gleichung (191) kénnen mit jedem 
Wert von a im allgemeinen drei Werte fir © = b, +b, bestimmt werden. Aus (190) ergibt sich 
dann S$ = b, bj. Aus © und $8 kann 6, und b, ermittelt werden, fiir die sich héchstens neun Werte 
ergeben. Der Wert ¢,/c, wird aus (192) bestimmt und damit ist auch (,/C. bekannt. Aus (193) bis 
(201) werden c, und ¢, errechnet, womit alle fiinf VektorgréBen bestimmt sind, die das sechs- 
gliedrige Kurbelgetriebe in seinen Abmessungen und in seiner Ausgangslage festlegen. 


6. SchluBbemerkung. Wie die Entwicklungen zeigen, kénnen mit Hilfe der Vektordarstel- 
lung in der Gaufschen Zahlenebene einfache Formeln fiir Winkelzuordnungen und Beziehungen 
ihrer Ableitungen aufgestellt werden, die zur Berechnung von Abmessungen an Kurbelgetrieben 
mit beliebiger Genauigkeit dienen kénnen und dabei das Zugrundelegen von Bedingungen zu- 
lassen, die bisher mit anderen Methoden nicht erfillt werden konnten. Die letzten beiden Ab- 
schnitte enthalten Entwicklungen, in denen die Anwendung der algebraischen Methode auf Pro- 
bleme der Getriebesynthese gezeigt wurde, die bisher als unlésbar galten. Die Darstellung ist 
einfach, wenn auch die Durchfiihrung der Berechnung ohne besondere Hilfsmittel nicht ohne 
weiteres méglich ist. Die einfachen Ergebnisse lassen jedoch erhoffen, daf sich weitere aufschlub- 
reiche Untersuchungen durchfiihren lassen werden. 


(Eingegangen am 31. Januar 1956) 
Anschrift des Verfassers: Dr.-Ing. K.-H. Sieker, Berlin-Schiachtensee, NiklasstraBe 62 
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Experimenteller Beitrag zur Entstehung turbulenter Strémung 
im Rohr 


Von J. Rotta 


1. Einleitung. Die turbulenten Stémungen gehéren heute zu den wesentlichsten Objekten der 
Strémungsforschung, da noch eine groBe Zahl offener Fragen zu beantworten ist. Als sehr be- 
deutend gelten die Fragen, die mit der Entstehung der Turbulenz zusammenhangen. Daher 
hat sich in neuerer Zeit das Interesse der Strémungsforschung viel den Gebieten zugewandt, 
in denen weder laminare noch voll turbulente Strémung herrscht. Za sehr aufschluBreichen 
Ergebnissen fiihrten experimentelle Untersuchungen der Strémungsstruktur bei der Platten- 
grenzschicht. 

Fur die Strémung im zylindrischen Rohr sind ahnliche Untersuchungen nicht bekannt 
geworden, obwohl diese Strémung von einfacherem Typus ist und manche grundsatzlichen 
Effekte dabei deutlicher hervortreten. Es ist die Tatsache, daB die Rohrstrémung nach einer 
gewissen Hinlauflange einem stationdaren Endzustand zustrebt, so da dann weiterhin die Lauf- 
lange als EinfluBgréBe fortfallt, die dieser Strémung Bedeutung als Studienobjekt zur Lésung 
grundsatzlicher Fragen verleiht. Dies war der AnlaB, Versuche in der Nahe der kritischen 
Reynoldszahl durchzufiihren, die hier mitgeteilt und diskutiert werden sollen. 


Eine grofe Menge von Erkenntnissen tiber das Auftreten turbulenter Str6mung im Rohr ist den 
Arbeiten L. Schillers! zu verdanken. Es ist bekannt, daB die GréBe der Reynoldszahlen, bei denen 
der Ubergang von laminarer zu turbulenter Str6mung im Rohr stattfindet, weitgehend durch die 
Strémungsverhaltnisse am Rohreinlauf bedingt ist. Aus der Versuchserfahrung weif man ferner 
seit langem, daB dieser Ubergang im allgemeinen nicht plétzlich vor sich geht, sondern daB es 
einen endlichen Bereich von Reynoldszahlen gibt, in denen die Strémung weder ganz laminar 
noch ganz turbulent ist. Aber uber die Ursachen dieses Verhaltens bestehen noch jetzt eine 
Anzahl von Unklarheiten. 


Gerade hinsichtlich der Rohrstrémung hatte die Tatsache, daB es zwei verschiedene Stré- 
mungsformen gibt, sehr zeitig das Interesse der Forscher angeregt. Vor allem hatte sich O. Rey- 
nolds? mit dem Problem beschaftigt und unter anderem festgestellt, daB man bei allmahlicher 
Erhéhung der Geschwindigkeit durch einen Bereich kommt, in welchem die Strémung abwech- 
wechselnd laminar und turbulent ist. Beim AusflieBen von Wasser aus einem Rohr wird dieses 
intermittierende Auftreten turbulenter Strémung, das fiir das Ubergangsgebiet typisch ist, 
durch Pendeln des frei ausflieBenden Strahles sichtbar. Aus diesen Pendelbewegungen wurde 
auf zeitliche Schwankungen der AusfluBgeschwindigkeit geschlossen. Die Erscheinung wurde 
von O. Tietjens® wie folgt interpretiert: 


,.Dabei entsprechen den Héchstwerten der Geschwindigkeiten immer diejenigen Augenblicke, 
in denen plétzlich innerhalb einer gewissen Rohrstrecke Turbulenz einsetzt. Durch den hiermit 
bedingten gréBeren Stro6mungswiderstand im Rohr nimmt die AusfluBgeschwindigkeit mit der 
Zeit ab, waihrenddessen die turbulent stroémende Fliissigkeitsmasse durch das Rohr hindurch- 
strémt. Sobald diese das Rohr zu verlassen beginnt, nimmt der Widerstand im Rohr wieder ab, 
was eine Zunahme der Geschwindigkeit zur Folge hat, bis-wieder teilweise Turbulenz einsetzt 
und der eben beschriebene Vorgang sich wiederholt.“‘ 


Unterstellt man die Richtigkeit dieser Deutung, so stellt sich die Frage, ob sich auch noch 
intermittierend turbulente Strémung zeigt, wenn die erwihnten Schwankungen der AusfluB- 
menge durch geeignete MaBnahmen verhindert werden, oder ob dann die Strémung einen anderen 
Charakter annimmt. Diese Fiage war der Ausgangspunkt der vorliegenden Versuche. Es war 
daher zuniachst eine Versuchsanlage zu schaffen, bei der die DurchfluBmenge vom Strémungs- 
widerstand des Versuchsrohres nicht beeinfluBt wird. 


L. Schiller, Stromung in Rohren, Handb. d. Exp. Physik, Bd. IV, 4. Leipzig 1932. 
0. Reynolds, Phil. Trans. Roy. Soc. London 1883 oder Scient. Papers Vol II, S. 51. 
O. Tietjens, Hydro- und Aeromechanix, Bd. JI, S. 38 ff. Berlin 1931. 
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2. Versuchseinrichtung und Durchfiihrung der Versuche. Mit Riicksicht auf beabsichtigte Hitz- 
drahtuntersuchungen sollte Luft als Strémungsmedium verwendet werden. Um die Grae der 
DurchfluBmenge von Schwankungen des DurchfluBwiderstandes unabhangig zu machen, war es 
naheliegend, die DurchfluBmenge vermittels einer Diise zu regeln, die mit Schallgeschwindigkeit 


durchstrémt wird. Ist Fp der Offnungsquerschnitt einer solchen Diise, so ist die in der Zeiteinheit 
durchstrémende Menge 
! 


: Denko SO are 
Mal) To VreePoes. ) 
Sie ist also auBer von Fp = x dj/4 nur von dem Druck Pp» der Dichte 9, der Luft vor der Diise 


und dem Verhaltnis der spezifischen Warmen x (x = 1,4 fiir Luft) abhangig. Ordnet man eine 
solche Diise vor dem Versuchsrohr an, so ist die DurchfluBmenge sicher von Druckschwankungen 
im Rohr unabhangig, solange der Druck hinter der Diise kleiner als der kritische Druck 


2 x—1 ; 
Prt = Pp | = 0,93 Pp (2) 
ist. 
Weiter sollte der Durchmesser des Versuchsrohres nicht zu klein gewahlt werden, um auch 


Geschwindigkeitsverteilungen messen zu kénnen. Um nun bei den kleinen Reynoldszahlen 
trotzdem Druckdifferenzen meBbarer GréBe zu erhalten, erwies es sich als zweckmabig, im Ver- 
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Abb. 1. Versuchsanlage. 


suchsrohr einen mdglichst niedrigen Druck p, zu haben. Es wurde deshalb die in Abb. 1 sche- 
matisch dargestellte Versuchsanordnung gewahlt, die an die im Institut vorhandene Unter- 
druckanlage (Pumpe und Unterdruckkessel) angeschlossen wurde. Die Luft wurde durch die 
erwahnte Diise, die bei Pos. VI angebracht war, aus dem Raum angesaugt und gelangte nach 
Durchstrémen der Beruhigungsstrecke mit eingebautem Gleichrichter in das eigentliche Ver- 
suchsrohr von insgesamt 10 m Lange und einem Innendurchmesser d = 30 mm, Der statische 
Druck p, in der Versuchsstrecke konnte durch mehr oder weniger starkes Offnen des Ventils bei 
Pos. I variiert werden. Es gelten folgende Beziehungen: 

Die Temperatur T, der Luft im Rohr sei1 gleich der Raumtemperatur Ty. Nach dem Gesetz 
der idealen Gase ist dann die Dichte der Luft im Rohr 


gree (3) 
| 01 Lo Po 


Zwischen durchstrémender Menge und der DurchfluBgeschwindigkeit U gilt die Beziehung 


. d? 
M=9,U7_. (4) 


1 Die Temperaturabminderung als Folge des Joule-Thompson-Effekts bei Durchstrémung der Diise ist 
vernachlassigbar. 
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Diese ist gleich der durch die Diise nach (1) gegebenen Menge. Die Rohr- Reynoldszahl ist 


jp) (5) 
bb 


Da nun die Zahigkeit ~ nur von der Temperatur der Luft abhangt, Q, U nach (4) ebenso wie M 
von p, unabhingig ist, so bleibt die Reynoldszahl Re unabhangig vom Druck p,, solange die 
Bedingung (2) erfiillt ist. Somit ist der Str6mungszustand in der Versuchsstrecke bei gegebener 
Diise und gegebenen Daten der Raumluft (pp, Qo, Ty) fiir jedes p, dbnlich. 

Der mit der Durchflu®geschwindigkeit U gebildete mittlere Staudruck q* errechnet sich zu 


Pier Msi wel Ay (6) 
tpg) 2 (d? 2/4)? 0, 2 (d? 2/4)? Qo Pi 


Dieser Staudruck und wegen der Ahnlichkeit der Strémung ebenso die mit anderen Geschwindig- 
keiten gebildeten Staudriicke und auch der Druckabfall im Rohr sind dem Druck p, umgekehrt 
proportional. Hier liegt einer der Vorteile dieser Versuchsanordnung, dafs die zu messenden 
Driicke durch Erniedrigung von p, erhéht werden kénnen, ohne daB sich dabei die Reynoldszahl 
andert. Es zeigte sich, daB bei einem praktisch realisierbarem Druck p, = 30mm Hg-Saule 
bei gleichen Abmessungen des Rohres und gleicher Reynoldszahl mit Luft als Strémungsmedium 
groBere Staudriicke erzielt werden konnten als es mit Wasser méglich gewesen ware. 


Im Normalfall sind der Druck Pp und die Dichte Op VOr der Diise gleich den entsprechenden | 


Werten der Luft im Raum. Aus den Gleichungen (1), (4) und (5) erkennt man dann, da fiir 
jede Reynoldszahl ein anderer Diisenéffnungsquerschnitt Fp erforderlich ist. Die Diisen bei 
Pos. VI waren deshalb leicht auswechselbar, und es war eine Anzahl verschiedener Groen vor- 
ratig. Eine noch feinere, stetige Anderung der Reynoldszahl konnte durch Beeinflussung des 
Druckes Pp Vor der Diise erméglicht werden. Dies wurde durch ein vor der Diise angebrachtes 
Drosselventil verwirklicht. Mit der nach dem Gasgesetz bei gleichbleibender Temperatur gelten- 
den Beziehung 


ee eD (7) 
Cn = Go pe 
geht (1) uber in 
1 
: 2 \ra One 
Mm = (735) Fy Pp Vrs (8) 


Die Messung des Druckes p, vor der Diise und der Raumtemperatur T, waren hiernach 


ausreichend, um die Reynoldszahlim Rohr zu bestimmen. Die zusatzliche Messung des statischen 
Druckes p, im Rohr erméglichte die Ermittlung der DurchfluBgeschwindigkeit U und des Stau- 
drucks q*. Diese Methode zur Bestimmung von Re und U ist sehr zuverlassig; Voraussetzung 
ist allerdings, daB die Versuchsanlage sorgfaltig abgedichtet ist, so daB nicht an Flanschverbin- 
dungen oder dgl. Nebenluft einstrémen kann. Die ausreichende Abdichtung bereitete keine 
Schwierigkeiten. Die benutzten Diisen wurden vorher durch eine Mengenmessung geeicht. 
Die GréBe von p, wurde durch ein Quecksilber-Barometer gemessen, ebenso der Druck der 
Raumluft. Zur Messung der Differenz p, — Pp zwischen den Driicken im Raum und vor der 
Diise wurde ein Prandtl-Manometer mit Quecksilberfillung benutzt. 


Zwischen der Beruhigungsstrecke und der Versuchsstrecke war ein 110 mm langer, allmah- — 


licher Ubergang vorgesehen. Das Kontraktionsverhaltnis war 1:11,1. An den StoBstellen des 
Rohres wurden die Innenkanten sehr sorgfaltig bearbeitet, so daB hier keine Stérungen der 
Strémung verursacht werden konnten. In diesem Zustand der Anlage blieb die Strémung bis 
zu Reynoldszahlen von Re = 18000 laminar. 

Um das Auftreten turbulenter Strémung méglichst nahe an die untere Grenze der Reynolds- 
zahlen zu bringen, wurden bei Pos. IV verschieden grobe Siebe oder Blenden eingebaut. Abb. 2 
zeigt die 3 Rohreinlaufe, bei denen die niedrigsten Reynoldszahlen fiir das Einsetzen der Turbu- 
lenz erzielt wurden. Bei der mit Nr. 1 bezeichneten Ausfiihrung war eine scharfkantige Lochblende 
hinter der Kontraktion eingebaut. Als nachster Schritt zur Verstarkung der Stérungen wurde 


die Beruhigungsstrecke fortgelassen und die Regeldiise direkt vor das Rohr gesetzt, Nr. 2. Auf | 


diese Weise wurde die Luft durch einen turbulenten Uberschallstrahl in das Versuchsrohr ein- 
gefiihrt. Bei Nr. 3 wurde die Luft durch 3 Uberschallstrahlen nahe der Rohrwand zugefihrt. 
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Mit dieser Versuchsanlage wurden Druckabfallmessungen, Hitzdrahtuntersuchungen und 
Staudruckmessungen durchgefiihrt. Die Hitzdrahtuntersuchungen dienten zum cecunen der 
Strémungsstruktur. Zur Verstarkung der von den Geschwindigkeitsanderungen a en 
Anderungen des elektrischen Stromes stand ein vorhandener Turbulenzverstarker a ae 
verstarker) nach H. Schuh! zur Verfiigung. Dieser wurde in Verbindung mit einem Elektronen- 
Oszillographen zur visuellen Beobachtung der Strémung benutzt. Da dieser Werecsthes Se 
schnelle Stromanderungen anspricht, war diese Methode zur Feststellung, ob laminare oder 
turbulente Strémung vorliegt, gut geeignet. Um aber auch die langsamer verlaufenden Ge- 
schwindigkeitsinderungen zu erfassen, wurde ein zwei- 
stufiger Gegentakt-Gleichspannungsverstarker gebaut, cpt ae be 
Abb. 3. Mit dessen Hilfe konnte der gesamte Geschwin- 
digkeitsverlauf iiber der Zeit von einem Schleifen-Oszillo- 
graphen sichtbar gemacht und aufgeschrieben werden, §// 


Fingang 1MQ 


an 


Abb. 2. 
Rohreinlaufe bei Pos. IV. Zweistufiger Gegentakt-Gleichspannungsverstarker. 


Bei den Staudruckmessungen und den Messungen des Druckabfalls wurden die Druck- 
differenzen mit einem Projektionsmanometer nach Betz gemessen. Daf dabei der statische 
Druck im Rohr sehr viel niedriger als im Raum war, hatte keinen unmittelbaren Einflu8 auf die 
Messung, wenigstens nicht, solange der Druck noch gréfer als der Dampfdruck des als MeB- 
fliissigkeit dienenden Wassers war. Zur Dimpfung der im Ubergangsgebiet zwischen laminarer 
und turbulenter Strémung auftretenden, heftigen Druckschwankungen wurden in die Leitungen 
zwischen MeS stellen und Manometer Kapillaren geschaltet, um eine zuverlassige Ablesung der 
zeitlichen Mittelwerte der Driicke zu erméglichen. Die Durchfiihrung genauer Druckmessungen 
war jedoch sehr miihselig und zeitraubend. Der hohe Widerstand der Druckleitungen hatte 
lange Einstellzeiten zur Folge und machte auBerdem die Anzeige selr empfindlich gegen minimale 
Undichtigkeiten in den Leitungen und im Manometer. In Verbindung mit den relativ grofen 
Luftvolumina unterschiedlicher GréBe iiber den beiden Fliissigkeitsspiegeln des Projektions- 
manometers ergab sich ferner eine groBe Empfindlichkeit gegen kleinste Anderungen des 


1 H. Schuh, Ein elektrisches Turbulenzmefgeraét, Deutsche Luftfahrtforschung UM 6607 (1944). 


262 Rotta: Experimenteller Beitrag zur Entstehung turbulenter Stromung im Rohr _Ingenieur-Archiv 


statischen Druckes p, (verursacht durch unregelmaBigen Lauf der Pumpe). Weitere Schwierig- 
keiten bereitete der Umstand, da8 in den Leitungen von Zeit zu Zeit Kondenswasser nieder- 
schlug, das die Kapillaren oder die Druckanbohrungen im Versuchsrohr verstopfte. Die Leitungen 
waren daher 6fter zu entliften. 

Die Staudruckmessungen zur Ermittlung der Geschwindigkeitsverteilungen, die mit einem 
iiblichen Pitotrohr durchgefiithrt wurden, muBten mit einer Korrektur fiir den EinfluB der Zahigkeit 
versehen werden, da die Reynoldszahl des Pitotrohres z. T. sehr klein war. Gemaf einer Arbeit 
von K. Wieghardt! wurde fiir diese Korrektur die Beziehung 


100 2) 


Tae a q (1 zu w dp (9) 


benutzt, wobei q der wirkliche Staudruck und q,,,,, der gemessene Wert des Staudrucks, u der 
zeitliche Mittelwert der Geschwindigkeit, dp der Durchmesser des Pitotrohres und y» die kine- 
matische Zahigkeit ist. Setzt man u? ~ 2 q/o, so ergibt sich fur den wirklichen Staudruck 


aoe 50 1? @ 
= IMess APH dj 


(10) 


Das Korrekturglied ist von der értlichen Geschwindigkeit unabhangig und damit fiir gegebenen 
Strémungszustand tibe: den Querschnitt des Rohres konstant. 


Die aus den Pitotrohrmessungen ermittelten Geschwindigkeitsverteilungen stimmten hin- 


sichtlich der DurchfluBmenge ; 
R 
M=2n0furdr (11) 
é 


mit der nach (8) sich ergebenden Menge mitunter nicht genau tiberein. Dies war insbesondere 


bei kleineren DurchfluBmengen der Fall. Da die Bestimmung von M nach (8) zuverlassiger er- 
scheint, sind die Ursachen fiir die Diskrepanz wahrscheinlich in den Unzulanglichkeiten der 
Pitotrohrmessungen = zu 
suchen. Neben zufalliger 
Streuung einzelner Mef- 
punkte sind wohl haupt- 
sachlich kleine unkontrol- 
lierbare Undichtigkeiten 
und die geringe, dem 
Beiz-Manometer _ eigene 
Schwankung des Null- 
punktes hierfiir verant- 
wortlich zu machen. Um 
diesen Einfliissen Rech- 
nung zu tragen, wurde 
eine weitere, kleine, tber 
den Querschnitt jeweils 
konstante Korrektur Aq 
angebracht, die so fest- 
gelegt wurde, daf die 
DurchfluBmenge nach 


Versuchsrohr 


laminar turbulent 5 aH 

stromende Menge strimende Menge (11) mit der nach (8) tiber- 
einstimmt. Die nach dieser 

Abb. 4, Getrenntes Auffangen turbulent und laminar ausstrémender Fliissigkeit und Zahlung Methode ermittelten Ge- 


der Strémungswechsel. ° : : : 
schwindigkeitsverteilun- 


gen bei laminarer Stromung brachten ausgezeichnete Ubereinstimmung mit der Theorie. 
Die Auswertungsmethode wurde daher als brauchbar erachtet und auf das Ubergangsgebiet 
der Str6mung angewandt. 

AuBer den Versuchen mit der beschriebenen Anlage wurden noch Versuche mit Glasrohren 
und Wasser als Strémungsmedium durchgefithrt, bei denen Beobachtungen des frei ausflieBenden 
Strahles vorgenommen wurden. Die Glasrohre hatten etwa 3 bis 6mm Innendurchmesser und 


‘ K. Wieghardt, Einige Bemerkungen iiber das Pitotrohr und die Zylindersonde, unveroffentlichter Bericht 
aus dem Kaiser-Wilhelm-Institut fiir Stromungsforschung (1946). 
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waren bis zu 1,50m lang. Das Wasser wurde einem 7m iiber dem Versuchsrohr gelegenen 
Vorratsbehalter entnommen. Der Strémungswiderstand des Versuchsrohres war klein gegeniiber 
dem dieser Fallhéhe entsprechenden Druck. Zwischen Rohr und Behiilter befand sich eae Danesel 
zur Regelung der DurchfluBmenge, deren Widerstand daher groB im Vergleich zum Rohrwider- 
stand war. So war auch bei diesen Versuchen die Gefahr ausgeschlossen, daB Schwankungen 
des Rohrwiderstandes Riickwirkung auf die DurchfluBmenge hatten. Bei einem Teil dieser 
Versuche wurde die Wurfweite des austretenden Strahles mit einem Theodoliten vermessen; 
die in der Zeiteinheit ausflieBende Wassermenge wurde mit MeBzylinder und Stoppuhr ermittelt. 


Die bekannten, schon erwahnten Pendelbewegungen des Strahles die im Bereich der kriti- 
schen Reynoldszahl auftreten, erméglichen es, die turbulent und laminar ausstr6mende Wasser- 
menge mit einem Doppeltrichter getrennt aufzufangen und zu messen, Abb. 4. Eine oberhalb 
der Trennwand im Trichter angebrachte Elektrode schlo& jeweils beim Ankommen der laminaren 
Str6mung einen Stromkreis zu einem Zahlwerk, so daB auch die Anzahl der Wechsel ermittelt 
werden konnte. Der Reiz dieser Versuchsmethode liegt in ihrer Einfachheit. Jedoch war einige 
Sorgfalt erforderlich, um reproduzierbare Versuchsergebnisse zu erhalten. Die Verteilung der 
Menge und die registrierte Zahl der Wechsel ist sehr von der Einstellung des Trichters zum Strahl 
abhangig. Es war notwendig, eine eigens hierfiir markierte Stelle des Trichters sehr genau zum 
Jaminaren Strahl auszurichten. AuffSerdem dndert sich in einem bestimmten Reynoldszahl- 
Bereich das Verhaltnis der Mengen so rasch mit der Reynoldszahl, daB ein MeBfehler in der 
Gesamtmenge die gleiche Abweichung von der mittleren Kurve ergibt wie ein 10mal so groBer 
Fehler des Verhaltnisses der getrennt aufgefangenen Mengen. Hine gewisse Streuung der Ver- 
suchsergebnisse muBte aus diesem Grunde in Kauf genommen werden. 
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Abb. 5. Druckabfallmessungen. 


3. Druckabfallmessungen. Die Untersuchung der Strémung im Bereich der kritischen Reynolds - 
zahl wurde mit einigen Messungen des Druckabfalls begonnen. Auf Abb. 5 sind die in der Versuchs- 
anlage nach Abb. 1 gewonnenen Ergebnisse mit den Rohreinlaufen nach Abb. 2 miteinander 
verglichen. Der dort iiber der Reynoldszahl Re aufgetragene Rohrwiderstandsbeiwert 


= PaP a ofl 
er = (12) 
oe 


wurde aus der gemessenen Differenz p, — p, der zeitlich gemittelten statischen Driicke an den 
Positionen II und I errechnet. Ax ist die Entfernung der MeBstellen voneinander. 


Bei Einlauftyp Nr. 1 findet der Anstieg von A infolge Turbulentwerdens der Strémung bei 
etwas gréBeren Reynoldszahlen statt als bei Einlauf Nr. 2, bei dem die Diise unmittelbar vor das 
Versuchsrohr gesetzt wurde. Die einstromende Luft fiihrt bei diesem Einlauf eine unvergleichlich 
héhere kinetische Energie mit sich als bei Einlauf Nr. 1. Bei Einlauf Nr. 3 hat die kinetische 
Energie der einstrémenden Luft etwa den gleichen Betrag, wie bei Kinlauf Nr. 2, nur wird jetzt 


die Luft mit ihrer Energie in Regionen eingefiihrt, die fiir die Turbulenzentstehung giinstiger 


sind. Dementsprechned beginnt der Anstieg von / bei noch kleineren Reynoldszahlen. 
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Diese Ergebnisse veranschaulichen zweierlei: 

1. In allen Fallen besteht trotz groBer Rohrlange ein endliches Intervall von Reynoldszahlen, 
in denen die Strémung nicht laminar, offenbar aber auch nicht voll turbulent ist. 

2. Die GréBe und Art der Stérung im Einlauf hat auch bei grobsten MaBnahmen noch eine 
Auswirkung auf die GréBe der Reynoldszahlen, bei denen sich der Wechsel von laminarer zu 
turbulenter Strémung vollzieht. Es bedarf wohl eines erheblichen Aufwandes, um die Strémung, 
die der unteren Grenze fiir diesen Wechsel entspricht, experimentell darzustellen, obwohl kein 
Zweifel besteht, daB es eine solche untere Grenze gibt. 


4. Hitzdrahtuntersuchungen. Im Bereich der kritischen Reynoldszahl angestellte Beobachtungen 
der Strémung am Ende des Rohres (bei Pos. I), bei denen die vom Hitzdraht gelieferten Strom- 
schwankungen auf dem Schirm eines Elektronen-Oszillographen sichtbar gemacht wurden, 
zeigten, daB das von AusfluBversuchen her bekannte intermittierende Auftreten turbulenter 
Strémung auch hier eintritt. Und zwar konnte die Erscheinung bei allen untersuchten Formen 
des Einlaufs festgestellt werden. Schwankungen der Durchflu&menge durch die Diise am Eintritt 


waren ausgeschlossen; es war aber noch zu priifen, ob nicht die Beruhigungsstrecke und das Rohr | 


ein Reservoir darstellte, das unter dem EinfluB der Kompressibilitat der Luft Ursache zu 
Schwankungen der DurchfluBmenge im Rohr gab. Dies ware denkbar, wenn die Schwankungen 
des Druckabfalls mit der GréBe des statischen Druckes in der Versuchsstrecke vergleichbar 


wiren. Bei den niedrigsten Werten des statischen Druckes p, = 30 mm Hg-S. betragen aber | 


die méglichen Schwankungen des Druckabfalls nur wenig mehr als 1% von p,. Die Hitzdraht- 
untersuchungen wurden jedoch meistens bei héheren Driicken (p, = 230 mm Hg-S.) ausgefuhrt, 
weil die Schwankungen dann wegen der kleineren Strémungsgeschwindigkeiten besser zu beob- 
achten waren. Hierbei waren aber die Schwankungen des Druckabfalls kleiner als 1°/o) des 


statischen Druckes. Damit steht fest, daB das Pendeln zwischen laminarem und turbulentem | 


Strémungszustand ohne Schwankungen der DurchfluBgeschwindigkeit einhergeht. Dies ist eine 
der wichtigsten Feststellungen der Untersuchungen. 


Abb. 6. Geschwindigkeitsschriebe in Rohrmitte bei Re = 2550. Rohreinlauf Nr. 1. Lauflange x/d = 322. Obere Kurve: Aufzeichnung nur 
rascher Geschwindigkeitsschwankungen (Hitzdrahtstrom durch Wechselstromverstarker geleitet); untere Kurve: Gesamtverlauf der Ge- 
schwindigkeit (Hitzdrahtstrom durch Gleichspannungsverstirker geleitet), Geschwindigkeit win m/sek.; Zeit t in sek. 


Die Eigentiimlichkeiten der Strémung werden durch eine Anzahl von Geschwindigkeits- 


schrieben veranschaulicht. Und zwar war die Versuchsanlage bei der Aufnahme dieser Schriebe, | 
ebenso wie bei allen weiteren in der Folge wiedergegebenen Versuchsergebnissen mit dem Rohr- | 


einlauf Nr. 1, also mit der Lochblende hinter der Kontraktion ausgeristet. 
Abb. 6 zeigt einen Schrieb! bei Re = 2550, bei dem sich der Hitzdraht in Rohrmitte bei 
Pos. I befand. Die Stromschwankungen wurden einmal iiber den Turbulenzverstarker geleitet 


(obere Kurve) und zum anderen iiber den Gleichspannungsverstarker (untere Kurve). Aus der | 


oberen Kurve sieht man deutlich das wechselweise Auftreten turbulenter und laminarer Strémung. 
Die untere Kurve gibt den tatsachlichen Geschwindigkeitsverlauf tiber der Zeit wieder und zeigt, 


da} jedesmal bei laminarer Strémung die Geschwindigkeit in Rohrmitte betrachtlich iiber dem | 


Niveau der turbulenten Strémung liegt. Wie man ferner sieht, ist die Dauer der einzelnen Peri- 
oden unregelmaBig; die bei laminarer Strémung auftretende Erhéhung der Geschwindigkeit ist, 
mit Ausnahme bei sehr kurzzeitigem Auftreten laminarer Strémung, jedoch immer die gleiche. 


Auf Abb. 7 sind entsprechende Geschwindigkeitsschriebe, unter Fortlassung der Tur- | 


bulenzverstaérker-Kurven, bei verschiedenen Reynoldszahlen dargestellt. Das unregelmaBige | 


Auftreten turbulenter Strémung wird mit abnehmender Reynoldszahl immer seltener. Bei Re = 


* Die Originalschriebe sind sehr lichtschwach und deshalb schlecht zu reproduzieren. Sie wurden daher 
photographisch etwa 4,5-fach vergréSert, dann mit der Hand nachgezogen und wieder verkleinert. 
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Pe. ist nur noch vereinzelt turbulente Strémung zu beobachten. Letzterer Schrieb erweckt den 
nschein, als ob das Wechseln jetzt von regelmaiBiger Dauer wire; dies ist jedoch nur ein Zufall 
Bei Reynoldszahlen kleiner als Re = 2200 ist die Strémung vollkommen laminar; die im Hinlaet 
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Abb. 7. Geschwindigkeitsschriebe. Gesamtverlauf der Geschwindigkeit in Rohrmitte bei verschiedenen Reyno(dszahlen, Robreinlauf Nr. 1. 
Lauflange x/d = 322. Geschwindigkeit ui n m/sek.; Zeit tin sek. 


7} 
7/R=02 


Abb. 8. Geschwindigkeitsschriebe. Gesamtverlauf der Geschwindigkeit bei Re =2550 in verschiedenen Abstanden von Rohrmitte. 
Rohreinlauf Nr. 1. Lauflinge x/d = 322, U = 4,27 m/sek. Geschwindigkeit win m/sek. ; Zeit tin sek. 


erzeugten Stérungen sind ganz abgeklungen. Bei Re = 2800 herrscht voll turbulente Strémung 
bei Pos. I. 

Wie das Strémungsbild an verschiedenen radialen Abstanden r/R von Rohrmitte aussieht, 
veranschaulicht Abb. 8 fiir Re = 2550. Aus diesen Geschwindigkeitsschrieben! muf man fol- 


1 Die Schriebe wurden nacheinander aufgenommen. Daher besteht keine zeitliche Zuordnung fiir dic 
turbulenten Partien bei verschiedenen r/R. 
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gern, daB sich bei der abwechsend laminaren und turbulenten Strémung auch, zumindest ange- 
nahert, die Geschwindigkeitsverteilungen fiir die jeweils zugehérige ausgebildete Rohrstrémung 
einstellt. Die zeitlich gemittelten Geschwindigkeitsverteilungen bei ausgebildeter Rohrstr6mung 
bezogen auf die DurchfluBgeschwindigkeit U, gibt Abb. 9 wieder. Diese Abbildung erlautert die 
auf Abb. 8 dargestellten Ergebnisse, namlich daB in Rohrmitte wahrend der laminaren Stromung 
eine héhere Geschwindigkeit als bei turbulenter herrscht. Weiter von der Mitte entfernt nehmen 
diese Unterschiede ab. Beir/R = 0,6 ist der Unterschied nur noch unbedeutend, und beir/R = 
0,8 ist die Geschwindigkeit bei laminarer Stromung merklich kleiner als bei turbulenter. 

Durch Auszihlen und Ausmessen der Schriebe 
wurde eine quantitative Analyse gewonnen, die 
allerdings keinen Anspruch auf groBe Genauigkeit 
erheben kann, da das Versuchsmaterial zu wenig 
umfangreich fiir eine statistische Auswertung ist. 
Zur quantitativen Beschreibung der Struktur der 
Str6mung wurden als statistische Merkmale ein 
Intermittenzfaktor y, der das durchschnittliche Ver- 
haltnis der turbulent strémenden Menge zur Ge- 
samtmenge ausdriickt, und eine Intermittenzzahl n, 
die die durchschnittliche Anzahl der Wechsel in der 
Zeiteinheit angibt, eingefiihrt. Unter einem Wechsel 
sei dabei eine volle Periode vom Beginn der lami- 
naren Strémung, einschlieBlich Ubergang zur tur- 
bulenten Strémung bis zum Wiedereinsetzen der 
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Abb. 9. Geschwindigkeitsverteilungen Abb. 10. Zur quantitativen Analyse 
bei ausgebildeter Rohrstrémung. der Hitzdrahtschriebe. 


laminaren Strémung verstanden. Ist N die Zahl der innerhalb der Zeit T beobachteten Wechsel, 
so ist die Intermittenzzahl n = N/T. Die Festsetzung der Grenze zwischen laminarer und 
turbulenter Strémung ist etwas willkiirlich und wurde entsprechend Abb. 10 vorgenommen. Der 
Intermittenzfaktor y ist dann definiert durch die Summe der Zeiten tr in denen turbulente 
Strémung herrrscht, dividiert durch die Beobachtungszeit T, also y = 3) t;/T. Die Inter- 
mittenzzahln, mit der DurchfluBgeschwindigkeit U und dem Rohrdurchmesser d dimensionslos 
gemacht, ist auf Abb. 11 tiber der Reynoldszahl aufgetragen. Darunter ist der Intermittenz- 
faktor y dargestellt Die dimensionslose Intermittenzzahl hat in der Nahe der Reynoldszahl 
Re = 2500 ihr Maximum, bei der y ~ 0,5 ist, d. h. bei der turbulente und laminare Strémung 
ungefahr gleich haufig sind. 

Bei den Hitzdrahtschrieben ist weiterhin noch das Verhalten an den Turbulenzfronten — 
wenn man die Ubergange von einem Strémungszustand in den anderen so bezeichnen will — sehr 
bemerkenswert. Der Ubergang vom turbulenten in den laminaren Zustand findet danach viel 


schroffer statt als der umgekehrte Ubergang vom Laminaren ins Turbulente. Fiir dieses Verhal- ” 


ten sind die Schriebe mit Hilfe des Turbulenzverstarkers, Abb. 6, sehr aufschluBreich. Unmittel- 
bar vor dem Beginn einer laminaren Partie sind die turbulenten Schwankungen der Geschwindig- 
keit besonders heftig, z. T. so heftig, da sie nicht mehr registriert wurden. Ganz plétzlich ver- 
schwinden sie dann. Die darauf folgende Anderung der Geschwindigkeit spielt sich innerhalb 
einer Zeitspanne ab, die wesentlich unter der Zeitkonstanten des Wechselstromverstarkers liegt. 
Der Verstarker reagiert daher mit einem Ausschlag, der dann gema der Zeitkonstanten ab- 
klingt. Der Geschwindigkeitswechsel bei dem Ubergang zu einer turbulenten Partie ist dagegen 
so langsam, da} die erwihnte Erscheinung nicht auftritt. Auch die Geschwindigkeitsschwan- 
kungen wachsen in diesem Fall von Null stetig an. 


3. Staudruckmessungen. Als Erganzung zu diesen Untersuchungen wurden Staudruckmessun- 
gen zur Ermittlung der zeitlichen Mittelwerte der Geschwindigkeitsverteilungen ausgefiihrt. Bei 
der Bestimmung der zeitlich gemittelten Geschwindigkeiten aus den gemessenen zeitlichen 
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Mittelwerten der Staudriicke ist der quadratische Mittelwert der Geschwindigkeitsschwankungen 
von Einflu8, der aber nicht gemessen wurde. Eine theoretische Abschatzung fiir die Lings- 
schwankungen beruht auf der von L. A. Sackmann! benutzten Annahme daB sich Partien mit 
voll ausgebildeter laminarer und voll ausgebildeter turbulenter Strémung miteinander abwech- 
seln und hinreichend scharf ineinander iiberwechseln. Die Grundlagen und Ergebnisse sollen hier 
eingefiigt werden, da sie zur Vertiefung der Vorstellung iiber das Wesen der intermittierend turbu- 
lenten Strémung beitragen. 
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Rohreinlauf Nr. 1. Lauflange x/d = 322. bei Re = 2500. 


Wahrend des Bruchteiles y der Zeit, in welcher turbulente Strémung herrscht, ist die drtliche 
Geschwindigkeit im Durchschnitt gleich derjenigen der ausgebildeten turbulenten Rohrstrémung, 
u = ur, und wahrend des Bruchteiles 1 — y, in welchem die Strémung laminar ist, nimmt wu den 
Wert der ausgegebildeten laminaren Rohrstr6mung an, u = uy. Damit ergibt sich der zeitliche 
Mittelwert der Geschwindigkeit zu 

u=yur+(l—y) uy. (13) 


Der quadratische Mittelwert der Geschwindigkeit ist dann 
. uw = yu, + (1—y) ud. (14) 


Bezeichnet man die Schwankungen der Geschwindigkeit mit Au = u — u, so folgt aus (13) und 
(14) fiir den quadratischen Mittelwert der Geschwindigkeitsschwankungen 


Au? = iw = (1—y) (uur). oc) 


teh, Ab Sacknann,..C. B. 224, (1948), S. 793, S.1326, C. R. 226 (1948) S. 1248, S. 1343, S. 1887, C. R. 227, 
(1948) S. 328, C. R. 239, (1954) S. 220, S. 389. 
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Mit den bekannten Geschwindigkeitsverteilungen bei ausgebildeter Rohrstrémung, Abb.9,und dem 
Intermittenzfaktor y nach Abb. 11 1aBt sich der quadratische Mittelwert fiir die Geschwindigkeits- 
schwankungen in Rohrachsenrichtung abschatzen. Auf Abb. 12 sind die sich so fur Rohrmitte 
ergebenden Werte iiber der Reynoldszahl aufgetragen und Abb. 13 gibt eine Vorstellung von der 
Verteilung iiber den Rohrhalbmesser bei Re = 2500. In die letztere Darstellung ist noch der 
Beitrag der turbulenten Schwankungen, der nach den Messungen von J. Laufer + abgeschatzt 
wurde, mit eingetragen. Die Schwankungen der Geschwindigkeit unterscheiden sich im vor- 
liegenden Fall sowohl nach GréBe als auch hinsichtlich der Verteilung iiber den Radius erheblich 
von denjenigen voll turbulenter Strémung. 
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Abb. 14. Zeitlicher Mittelwert der Geschwindigkeitsverteilung bei intermittierendem Auftreten turbulenter Strémung. 
Rohreinlauf Nr.1. Lauflange x/d =322. 


Da diese Schwankungen im allgemeinen wesentlich gréBer als diejenigen voll turbulenter 
Rohrstrémung sind, ist es erforderlich, bei der Auswertung der Staudruckmessungen wenigstens 
die Langsschwankungen zu beriicksichtigen. Der zeitliche Mittelwert der Geschwindigkeit er- 
rechnet sich dann aus dem gemessenen und korrigierten zeitlichen Mittelwert des Staudrucks 


naherungsweise zu 
— 2 [aaa 
b= 0 q Tie ta) (16) 


Auf Abb. 14 sind die so an der MeBstelle Pos. I ermittelten Geschwindigkeitsverteilungen fiir 
verschiedene Reynoldszahlen dargestellt. Man erkennt den allmahlichen Ubergang vom Profil der 
turbulenten Strémung in das Parabelprofil der laminaren Strémung, wenn die Reynoldszahl ab- 
nimmt. Die Geschwindigkeitsverteilungen an den Mefstellen Pos. II und III sind dbnlich. 
Nur erstreckt sich der Ubergang von turbulenter zu laminarer Strémung an den Stellen naher 
dem Einlauf ber ein gr6Beres Intervall von Reynoldszahlen. Dies veranschaulicht Abb. 15, auf 
der iiber der Reynoldszahl die Geschwindigkeit in Rohrmitte aufgetragen ist, wie sie sich se den 
Staudruckmessungen an den Mefstellen Pos. I, II und III ergibt. Es ist kein Anzeichen dafiir 
opee is da bei Pos. I ein stationarer, d. h. von der Rohrlange unabhangiger Endzustand 
erreicht ist. 


Die Staudruckmessungen erganzen also das aus den Hitzdrahtuntersuchungen gewonnene 
Bild durch das Ergebnis, daf die statistischen Higenschaften der Strémung sich in der Nahe der 
kritischen Reynoldszahl mit der Lauflange andern und selbst nach einer Entfernung von mehr als 
300 Durchmessern vom Einlauf noch keinen Endwert erreicht haben. 


1 J. Laufer, The structure of turbulence in fully developed pipe flow, NACA T. R. No. 1174 (1954) 
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6. Untersuchungen zur Wurfweite eines ausflieBenden Wasserstrahles. Bei den bisher bekannt 
gewordenen Untersuchungen tiber das intermittierende Auftreten turbulenter Strémung lieB man 
Wasser aus einem Behalter durch ein Rohr strémen. Aus den beobachteten Pendelbewegungen 
des ausflieBenden Strahles wurde auf zeitliche Schwankungen der DurchfluBgeschwindigkeit U 
geschlossen, wie schon erwahnt wurde. Von dieser Voraussetzung ausgehend, hat in neuerer Zeit 
L. A. Sackmann! die Rohrstrémung im Bereich der kritischen Reynoldszahl untersucht und mee 
lysiert. Wenn man von den Folgerungen betreffs der Schwankungen von U absieht, so decken 
sich die von ihm mitgeteilten Beobachtungen mit denen der hier beschriebenen Versuche weit- 
gehend. So wurde festgestellt, daB das Pendeln des AusfluGstrahles im wesentlichen zwischen 
zwei Extremlagen stattfindet und in unregelmaBiger Folge vor sich geht. Auch wurde von Sack- 
mann beobachtet, daB der Ubergang von turbulenter in laminare Strémung viel rascher erfolgt 
als der umgekehrte Wechsel von laminarer in turbulente Strémung. Es ist daher zu vermuten, daB 
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Abb. 15. Zeitlicher Mittelwert der Geschwindigkeit in Rohrmitte bei intermittierendem Auftreten turbulenter Strémung. Rohreinlauf Nr. 1. 


die Strémung dort von der gleichen Art wie bei den hier beschriebenen Versuchen mit Luft war. 
Diese Vermutung wird noch dadurch gestiitzt, daB nach den Hitzdrahtschrieben die turbulenten 
und laminaren Partien so haufig miteinander abwechseln, da} bei der von Sackmann benutzten 
Rohrlange niemals abwechselnd laminare und turbulente Strémung uber die ganze Rohrlange 
hergestellt wird. Im Bereich der groBten Intermittenzzahl befinden sich bei einer Rohrlange von 
200 Durchmessern etwa 6 bis 7 turbulente Partien gleichzeitig im Rohr. Es ist undenkbar, da die 
erzeugten Schwankungen des Rohrwiderstandes iiberhaupt nennenswerte Schwankungen der 
DurchfluBmenge zur Felge haben kénnen. Vor allem sind aber die unregelmafige Folge der 
Wechsel und die RegelmafBigkeit der Pendelamplituden unter solcher Voraussetzung nicht mit- 
einander vereinbar . Dies gibt Veranlassung, die Wurfweite eines ausflieBenden Strahles naher zu 
untersuchen. 

Man kann aus der Wurfweite des Strahles einer aus einem Rohr flieBenden Flissigkeit durch- 
aus auf die DurchfluBgeschwindigkeit U im Rohr schlieBen; es sind aber mehrere Nebeneinfliisse 
dabei zu beriicksichtigen. Insbesondere ist die sich unter der Wirkung der Reibung an der Rohr- 
wand einstellende Geschwindigkeitsverteilung tiber den Querschnitt von Bedeutung. Nach dem 
Verlassen des Rohres fallt die Wandreibung fort, und die Flissigkeit hat unter der Wirkung von 
Zahigkeit bezw. Turbulenz und Kohasion die Tendenz sich gleichmafig tiber den Querschnitt zu 
verteilen. Beim Ubergang zu dieser gleichmaBigen Geschwindigkeitsverteilung bleibt die in der 


Zeiteinheit austretende Menge M und der durch die Flache A der Kontrollflache, Abb. 16, tre- 


tende Impuls = 
D2. [adr (17) 


1 siehe FuBnote 1 von S. 267. 
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erhalten. Aus dem Impuls durch die Flache B 
J =», M (18) 


ergibt sich die fiir die Wurfweite des Strahles maBgebende Geschwindigkeit v, zu 


wi NP if 
0 


Fiir das Parabelprofil der laminaren Rohrstrémung ergibt sich hiernach v, = = U, wahrend bei 
turbulenter Strémung v, praktisch gleich U ist (v,; ~ 1,02 U). 

Bei den Abmessungen und Ge- 
schwindigkeiten, wie sie bei den in 
Frage stehenden Versuchen mit Wasser 
vorkommen, spielen auch die Kapillar- 
spannungen eine nicht zu ubersehene 
Rolle. Fiir eine Uberschlagsrechnung 
J wird einfach die an der Rohrmiindung 
Abb. 16. Anderung der ee ee ae phe beim ausflieBenden auf den Strahl ausgeiibte Kraft K von 

eon dem Impuls nach (18) abgezogen. Dies 
ist streng nur fiir sehr grofe Fallhéhe h 
zulassig. Zur Abschatzung von K selbst 
wird angenommen, da® der Durch- 
messer des Strahles an der Miindung 
gleich dem Innendurchmesser d des 
Rohres ist und sich mit dem Abstand 

Worfweite: «4 ye Messungen von der Miindung nicht andert. Die 
tf es) hid tam. turb. Dia Kapillarspannungen in axialer Rich- 
i | Tie PAE RN tung geben die Kraft Cad, wobei C 
os ia ls die Kapillarkonstante ist. Durch die 
ye1TNU ie is : : . 

Kapillarspannungen in Umfangsrich- 
tung steht die Fliissigkeit im Strahl 
i unter dem Uberdruck p; = 2 C/d: die 
SSS qa sich hierdurch ergebende Kraft pia d?/4 

s ™ muf yon der Kraft der axialen Span- 
nungen abgezogen werden. Somit ist 


die Kapillarkraft an der Miindung 


B 


turbulent 
08 


| | Cad 
Ve == K=-=, (20) 


\ Die einfache, aus dem Fallgesetz 
\ folgende Beziehung fiir die Wurfweite 
des Strahles lat sich nun benutzen, 
wenn die genannten Einfliisse durch 
einen Faktor k beriicksichtigt werden. 
(yay Fir die Wurfweite s ergibt sich dann 


OW 


O2 


0 007-002 ~—S—Ci08.—=—=<(‘«iC SC*<C«i OG SC«iSC<C*‘«iC SS Th 
Saul (21) 
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sek cm 
Atb. 17. Wurfweite des ausflieBenden Wasserstrahles. 

wobei h die Fallhéhe und g die Erd- 

beschleunigung ist. Nach den gemachten Abschatzungen ergibt sich fiir den Faktor k 


aC 
Ue lds 


(22) 


Die GréBe k wird weiter beeinfluBt durch die Lauflange, die zur Umwandlung in die gleichmabige 
Geschwindigkeitsverteilung nach dem Austritt erforderlich ist, durch die Wirkung der Erd- 


anziehung auf die Geschwindigkeitsverteilung kurz vor dem Austritt und durch die Reibung des 
Strahles mit der umgebenden Luft. 
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Zur Nachpriifung der Uberlegungen wurden einige Versuche durchgefiihrt. Die hieraus er- 
rechneten Werte des Faktors k sind auf Abb. 17 iiber (U ydy ‘aufgetragen. Zum Vergleich wurden 
die theoretischen Werte nach (22) mit eingezeichnet, denen die Kapillarkonstante C — 74 
< 10~* g/cm zugrunde liegt. Die MeBergebnisse bei turbulenter Strémung entsprechen sehr gut der 
Theorie; dagegen ergaben sich bei laminarer Strémung kleinere k-Werte als nach (22). Die we- 
sentliche Ursache ist, daB die Annahmen fiir die Abschatzung der Kapillarkrafte an der Miindung 
zu grob sind. Abb. 18 zeigt die Form des Strahles an der Miindung. Die Kapillarspannungen 
greifen am AuBendurchmesser D an; unter dem Einflu® der Erdbeschleunigung tritt zudem eine 
unsymmetrische Verformung auf. Fiir den vorliegenden Zusammenhang ist aber das Wichtigste 
die tatsachlich beobachtete 
Zunahme von k beim Uber- 
gang vonturbulenter zu lami- 
narer Strémung. 


Aus diesem Ergebnis ist zu 
folgern, daB die Pendelbewe- 
gungen des Ausflu®Bstrahles 
bei intermittierend auftreten- 
der, turbulenterRohrstré6mung 
keineswegs unbedingt zeitliche 
Schwankungen der Durchflub- 
geschwindigkeit U bedeuten 
mussen. Unter der Annahme, 
daB bei den Versuchen von 
Sackmann keine Schwankun- 
gen vonU auftraten, entspricht 
das Verhaltnis der dort er- 
mittelten Geschwindigkeiten 
der turbulenten bzw. lami- 
naren Partien zur Durchflub- 
geschwindigkeit dem Faktork 
in (21). Die so aus der Arbeit 
Sackmanns ermittelten k-Wer- 
te wurden in Abb. 17 mit ein- 
getragen und entsprechen im 
wesentlichen den ubrigen Er- 
gebnissen. Hiernach scheint 
Ess dal auch bei den Versuchen Abb. 18. Form des ausflieBenden Wasserstrahles an der Mitndung, VergréBerung 2,5:1. 
von Sackmann keine Schwan- a) Turbulenter Strahl (U |/d)~! = 0,012 (sekem=3/2); 
kungen der DurchfluBmenge b) Laminarer Strahl (U \/d)~! = 0,015 (sekem-®/2); 
auftraten. Die Pendelungen c) Laminarer Strahl (Ud)! = 0,052 (sekem=#/2). 
des Strahles werden durch die 
unterschiedliche GréBe des Impulses der Geschwindigkeitsverteilungen verursacht. Mit dieser 
Feststellung ist dann der Widerspruch mit den hier mitgeteilten Versuchen aufgeklart. 


a 


7. AusfluBversuche. Die Staudruckmessungen gaben den Hinweis, dafs nach einer Lauflange 
von einigen hundert Rohrdurchmessern noch kein stationarer Endzustand erreicht ist. Dies legt 
die Vermutung nahe, da die Anderung des Strémungszustandes mit der Rohrlange ein sehr 
wesentliches Merkmal ist, das fiir die Mechanik dieser Strémung vielleicht entscheidend ist. Da- 
her schien es wiinschenswert, zur Klarung dieser Frage weiteres Versuchsmaterial zu beschaffen. 


Nachdem im vorigen Abschnitt klargestellt worden war, daB die Pendelungen des aus einem 
Rohr ausflieBenden Wasserstrahles durch die unterschiedlichen Geschwindigkeitsverteilungen bei 
turbulenter und Jaminarer Str6mung verursacht werden, bietet sich die Méglichkeit, durch ge- 
trenntes Auffangen der turbulent und laminar ausstrémenden Wassermenge und durch Zahlen 
der Pendelbewegungen die fiir die Kennzeichnung der Struktur der Strémung wesentlichen 
GréBen des Intermittenzfaktors y und der Intermittenzzah! n auf verhaltnismafig einfache Weise 
zu ermitteln. Das anfanglich etwa 1,50 m lange Versuchsrohr von ungefahr 3 mm Innendurch- 
messer wurde nach jeder Versuchsreihe um ein Stiick verkiirzt. Auf diese Weise wurde ein Uber- 
blick iiber die Veranderung von und n mit der Lauflange x der Strémung erhalten. Die Ergeb- 
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nisse sind auf Abb. 19 und 20 iiber der Reynoldszahl dargestellt. Die Ordinaten-Nullpunkte wur- 
den fiir die verschiedenen x/d-Werte verschoben, um ein Uberdecken der MeBpunkte der ver- 
schiedenen Versuchsreihen zu vermeiden. Durch die Mefpunkte wurden geglattete Kurven ge- 
legt. Die Ergebnisse der drei zuerst durchgefiihrten MeBreihen (fiir x/d = 483, 435 und 387) sind 
hinsichtlich der Intermittenzzahl teilweise fehlerhaft!, Der mutmafliche Verlauf in dem frag- 
lichen Gebiet wurde in Abb. 20 und ebenso in den weiteren hieraus hergeleiteten Abbildungen 

. gestrichelt gekennzeichnet. 
Die geglatteten Kurven der Abb. 19 und 20 sind unter Vermeidung der Koordinaten Null- 
- punktsverschiebung noch einmal in Abb. 21 und 22 aufgetragen. Das Ergebnis zeigt, dab iiber 
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Abb. 21. Intermittenzfaktor als Funktion der Reynoldszahi! und der Lauflange. Gegliattete Kurven von Abb, 19. 
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Abb. 22. Dimensionslose Intermittenzzahl als Funktion der Reynoldszahl und der Lauflinge. Geglattete Kurven von Abb. 20. 


den ganzen Bereich der untersuchten Rohrlangen bis zu 500 Durchmessern eine stetige Anderung 
der Strémungsstruktur festzustellen ist. Der Bereich von Reynoldszahlen, in welchem inter- 
mittierende Strémung auftritt, wird mit der Rohrlange immer schmaler und verschiebt sich zu 
kleineren Werten von Re. Die Intermittenzzahl] nimmt mit der Rohrlainge ab. 

Quantitativ stimmen die Ergebnisse nicht genau mit der Auswertungen der Hitzdraht- 
schriebe iiberein, obwohl am Einlauf auch eine Ringblende gemaB Nr. 1 von Abb. 2 angebracht 


1 Infolge unzulanglicher Lagerung des Versuchsrohres waren Schwingungen moglich, die ein vorzeitiges 
Zerfallen des laminaren Strahles unter der Wirkung der Kapillarkrafte zur Folge hatten. Die Zahleinrichtung 
reagierte aber nicht nur auf die Wechsel zwischen laminarer und turbulenter Strémung, sondern sprach auch 
auf einzelne Tropfen an. Daher ist die gemessene Intermittenzzahl bei diesen Reihen zu gro8, allerdings nur 
in den Bereichen, in denen laminare Partien von nicht zu kleiner Lange auftraten. Die Ursache wurde zu spat 
erkannt; eine Wiederholung der Versuche war nicht méglich. 
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war. Die Bedingungen der Zustrémung zur Blende sind alleidings in beiden Fallen verschieden; 
darin mag ein Grund fiir die unterschiedlichen Ergebnisse liegen. Ein weiterer Grund kann in der 
Methode liegen: Die Abgrenzung zwischen turbulenten und laminaren Parten ist bei beiden 
Methoden bis zu einem gewissen Grade willkiirlich festgelegt worden. Hinsichtlich der grund- 
sitzlichen SchluBfolgerungen sind diese Unterschiede jedoch ohne Bedeutung. 


8. Uberlegungen zu den Ubergiingen an den Turbulenzfronten. Aus der Tatsache, daB in der 
Versuchsanlage, Abb. 1, mit glattem Rohreinlauf bis zu sehr groBen Reynoldszahlen laminare 
Strémung besteht, ist zu schlieBen, daB die bei kleineren Reynoldszahlen beobachtete Turbulenz 
ihren Ursprung ausschlieBlich im gestérten Rohreinlauf hat. Unter dieser Voraussetzung kann 
man das intermittierende Auftreten turbulenter Strémung nur so deuten, daB Partien laminarer 
und turbulenter Strémung durch das Rohr hintereinander herflieBen. Da mit den angewandten 
Versuchsmethoden keine direkten Anhaltspunkte fiir die Fortbewegungsgeschwindigkeit der 
turbulenten Partien gewonnen werden konnten, sollen die Strémungsverhdltsnisse an den Turbu- 
lenzfronten einer theoretischen Betrachtung unterzogen werden, um so aus den Versuchsbeob- 
achtungen weitere Schlisse zu ziehen. 

Zur Klarstellung der Begriffe sei die Strémung von einem Beobachter betrachtet, der sich mit 
einer turbulenten Partie fortbewegt. Ihm stellt sich die Abgrenzung zur vorauslaufenden, lami- 
naren Partie als ,,vordere Turbulenzfront” dar; dieser Ubergang ist ein allmahlicher, wie aus den 
Aufzeichnungen der Geschwindigkeiten hervorgeht. Die Begrenzung zur nachfolgenden, lami- 


naren Partie erscheint dagegen als ,,hintere Turbulenzfront“ und ist ziemlich scharf begrenzt. 

Dem sich mit der Lauf- 
geschwindigkeit V der Turbulenz- 
front bewegenden Beobachter er- 
scheint die Strémung in der Um- 
gebung der Front in zeitlichen 
Mittelwerten stationar, wenn V 
konstant ist. Deshalb werde vor- 
ubergehend ein Koordinaten-Sy- 
stem eingefiihrt, in welchem sich 
das Rohr mit V nach rechts be- 
wegt und die Turbulenzfront ruht, 
Kontroliflache Abb. 23. Die Turbulenzfront 
werde von einer zylindrischen 
Kontrollflache eingehillt, deren 
Mantel dicht an die Rohrwandung grenzt, und deren Grenzflachen so weit von der Turbulenz- 
front entfernt sind, das dort ausgebildete laminare bezw. turbulente Strémung herrscht. Die 
GréBen am Ort der Grundflache auf der laminaren Seite werden mit dem Index L, diejenigen 
auf der turbulenten Seite mit T gekennzeichnet. Die Strémungsgeschwindigkeit in diesem 
System u =u sei—V. Zur Unterscheidung der Laufgeschwindigkeiten der vorderen und 
hinteren Turbulenzfront werden die Indizes 1 und 2 eingefiihrt. 

An der hinteren Turbulenzfront wirken auf die Kontrollflache folgende Krafte: An den Grund- . 
flachen greift der statische Druck p, der iiber die Grundflache konstant vorausgesetzt wird, und 
der Impuls g u* der Strémung an; am Mantel wirkt die Wandschubspannung 1). Damit ergibt 
sich das Kraftegleichgewicht in x-Richtung 
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Abb, 23. Zum Ubergang an den Turbulenzfronten. 
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wR? (pr— pr) +22 ef (uy — wy) r dr — 2 at Rf x dx =0. (23) 
Da nun . 
R R 
2x fuprdr = 20 f upr dr =a RU (24) 
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ist, folgt aus (23) nach Wiedereinfiihrung von w statt w fiir das Impulsgleichgewicht die von der 
Laufgeschwindigkeit V, unabhangige Beziehung 
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Bei Anwendung auf die vordere Turbulenzfront sind lediglich die Indizes L und T zu vertauschen. 


. Wandschubspannung stets positiv. Auch 
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1 
Da in der laminaren Strémung das Impulsintegral 2 | (a) a d (x) groBer (= 1,333) als bei tur- 
0 


bulenter Strémung ist (~ 1,02), so entsteht beim Ubergang jeweils eine Druckaénderung. Vom 
Standpunkt des ruhenden Beobachters zeigt sich beim Ankommen einer turbulenten Partie 
eine Druckerhéhung, bei Ankunft einer laminaren Partie eine Drucksenkung. 

Durch die beim Impulsgleichgewicht beriicksichtigten auBeren Krifte wird dem von der 
Kontrollflache umschlossenen Volumen mechanische Energie zugefiihrt, die innerhalb dieses 
Volumens entweder in Warme (Dissipation) oder in kinetische Energie turbulenter Geschwindig- 
keitsschwankungen umgewandelt wird. Diese in der Zeiteinheit umgewandelten Energien werden 
im folgenden mit Ey (Dissipation) und E, (Schwankungen) beziehungsweise bezeichnet. Ent- 
sprechend Gl. (23) gilt jetzt 


R oor : c 
a FR? (p, —p,)(U—V,) +20 | (ui — uy) rdr + 2a RV, [to dx—E,—E,=0. (26) 
0 xr 

Dies 148t sich mit (24) und Wiedereinfiihrung von u umformen, wobei noch die Differenz der 
Driicke pr — pr mit (25) eliminiert werden kann. Es ergibt sich dann eine Gleichung von der 
Form 

: Sire ; ; 

E, +22 RU{1t)dx —E,—E, =0. (27) 

. *L 

Das erste Glied E, stellt den Energiebetrag dar, der dem von der Kontrollflache umschlossenen 
Volumen durch Anderung der Geschwindigkeitsverteilung zu- (oder ab-) gefiihrt wird. Dieser 
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und ist von der Laufgeschwindigkeit V, der hinteren Turbulenzfront abhangig. Er laft sich aus 
den bekannten Geschwindigkeitsverteilungen fiir ausgebildete turbulente und laminare Rohr- 
strémung Abb. 9, errechnen und ist auf 


Abb. 24 dimensionslos tiber dem Ge- 
schwindigkeitsverhaltnis V,/U aufgetra- 


gen. E, nimmt ab, wenn V,/U steigt 
und wird negativ, wenn V, die Grobe 
1,05 U ubersteigt. 

Wahrend E; sowohl positiv als auch 
negativ sein kann, ist die durch das zweite 
Glied in (27) reprasentierte Arbeit der 


Ey und E, sind immer positiv, denn so- 
wohl die Umwandlung von mechanischer 
Energie durch Reibung in Warme als 
auch die Erzeugung von Energie un- 
geordneter Turbulenzbewegung aus der 
geordneten, laminaren Strémung sind Abb. 24. Zur Energiebilanz an der hinteren Turbulenzfront. 

nicht umkehrbar. Der in turbulente Bae Energiegewinn durch Umformung der zeitlich gemittelten Geschwindig- 
Geschwindigkeitsschwankungen verwan- _,_‘**itsverteilungs ae 

alte Energiebetrag ae E, =in Turbulenzschwankungen umgesetzte Energie. 


’ R 
Ei, = 7 0 f (u’2 + v2 + w’?) (up — Vy) rdr. (29) 
0 


Eine Abschatzung dieser GréBe ist unter Zugrundelegung der bei ausgebildeter Rohrstromung 
auftretenden Schwankungen mdglich, die allerdings bei wesentlich gréBeren Reynoldszahlen 
gemessen wurden!. Diese Gréfe ist auf Abb. 24 ebenfalls dargestellt. 


1 Siehe FuBnote 1 von S. 268. 
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Die Betrage der Wandschubspannungsarbeit und der Energiedissipation hangen von der | 
Lange xp — xz, des Kontrollzylinders ab. Da bei ausgebildeter Rohrstrémung, ob laminar oder | 
turbulent, die beiden Betrage einander gleich sind, so hangt beim vorliegenden Problem ihre | 
Differenz von der Linge der Strecke ab, innerhalb der sich die Umwandlung vollzieht. Ist diese | 
Lange gleich Null, d. h. findet der Ubergang plotzlich statt, so kann man sich die Kontrollflache | 


} 


ganz dicht an die Turbulenzfront herangeriickt denken, so daB die Betrage der Wandschub- | 


spannungsarbeit und der Dissipation vernachlassigbar klein werden. In diesem Fall miissen sich 


die durch Umwandlung der Geschwindigkeitsverteilung gewonnene Energie E,, und die in kine- | 


tische Energie turbulenter Schwankungen verwandelte Energie E, das Gleichgewicht halten. 


Auf Grund dieser Uberlegungen ist ein stoBartiger Ubergang an der riickwartigen Turbulenz- | 
front vom Standpunkt des Energiehaushaltes miglich, wenn die Laufgeschwindigkeit V, der | 
Front kleiner oder zumindest nur unbedeutend gréBer als die DurchfluBgeschwindigkeit U ist. 
Die Versuche zeigen, daB die Geschwindigkeitsverteilung sich wohl rasch, aber doch nicht | 


schlagartig andert. Auf der Strecke, in der die Umwandlung erfolgt, wird sicher die Dissipation 


gréBer als die Arbeit der Wandreibung sein, so da hier noch ein Teil kinetischer Energie ver- | 
braucht ‘wird. Trotzdem ist es nicht wahrscheinlich daB die Geschwindigkeit V, dieser Turbu- | 


lenzfront sehr viel von der DurchfluBgeschwindigkeit abweicht; ganz unwahrscheinlich ist es 
aber, daB V, wesentlich gréBer als U ist. 


Fiir die Verhaltnisse an der vorderen Turbulenzfront sind in (28) wieder die Indizes L und T | 
zu vertauschen und V, fiir V, einzufiihren. E,, wechselt also einfach das Vorzeichen. Wenn V, 


kleiner oder nur unwesentlich gréBer als U ist, werden E, negativ und E,=0. Der zur Um- | 


wandlung der Geschwindigkeitsverteilung erforderliche Antrieb kann dann nur aus der Arbeit 


der Wandreibung kommen. Das geht so vor sich, daB in x-Richtung die Geschwindigkeits- | 


schwankungen abnehmen, die Verteilung der zeitlich gemittelten Geschwindigkeit sich aber 
langsamer andert. Damit gewinnt dann die Wandreibungsarbeit das Ubergewicht iiber die 
Dissipation. So findet in x-Richtung eine stetige Umwandlung der Geschwindigkeitsverteilung 


statt, die energiemaBig keinesfalls auf einer infinitesimalen Lange des Rohres méglich ist. Der | 
Vorgang ist offensichtlich dem der laminaren Einlaufstrémung im Rohr verwandt. Dies be- | 


statigen die Geschwindigkeitsschriebe sehr deutlich. 

Dieser Ubergang 1a®t sich unter vereinfachenden Rechnungs-Annahmen auch quantitativ 
verfolgen. Es wird dabei vorausgesetzt, die Umwandlung der Geschwindigkeitsverteilung werde 
durch plétzliches Verschwinden der turbulenten Scheinspannung — @ u’ v’ ausgelist. Ferner 
wird der Druck p an jeder Stelle x tiber den Querschnitt konstant angenommen, so daf sich die 
Impulsgleichung (25) in differentiierter Form schreiben laBt 


2 | (Bh a4] -— Be. 09 


0 Q ar) : (31) 
AuBerdem wird fiir die Verhaltnisse in Rohrmitte noch die Navier-Stokessche Bewegungsgleichung | 
OU», d 0? 
r—0 


hinzugenommen, mit deren Hilfe der Druckgradient dp/dx in (30) eliminiert werden kann. Die 
Integration der so entstehenden, gewéhnlichen Differentialgleichung erster Ordnung liefert die 
Lésung des Problems. Fiir die Geschwindigkeitsverteilung wird angenommen, da sie sich im 
Ubergangsgebiet aus den Verteilungen der ausgebildeten laminaren und turbulenten Stré6mung 
linear superponieren laBt: 

u =a uz + (1l—a) ur. (33) 
Die Ergebnisse dieser Rechnungen wurden umgerechnet auf das rohrfeste Koordinatensystem. 
Sie reprasentieren dann einen zeitlich verlaufenden Vorgang (s. Abb. 25), und sind somit mit 
den Hitzdrahtschrieben unmittelbar vergleichbar. Die groBe Ahnlichkeit aire Darstellung mit 
den Geschwindigkeitsschrieben spricht fiir die Richtigkeit der angestellten Uhcaseyi eee Der 
Abb. 25 ist noch das Ergebnis eines typischen Geschwindigkeitsschriebes beigegeben. Diesem 
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Vergleich nach zu urteilen, dirfte die Laufgeschwindigkeit V, der vorderen Turbulenzfront sicher 
auch nicht viel von U abweichen. 

: Die hier gezogenen Schliisse gipfeln darin, daB die experimentell beobachtete Form der Uber- 
ginge an den Turbulenzfronten vom Standpunkt des Energiehaushaltes Laufgeschwindig- 
keiten der Turbulenzfronten bedingen, die mit der DurchfluBgeschwindigkeit U der Gra fen- 
ordnung nach iibereinstimmen. Das an einem festen Ort des Rohres beobachtete inter- 

- mittierende Auftreten turbulenter Strémung entsteht dadurch, daB abwechselnd Partien tur- 
- bulenter und laminarer Str6mung durch das Rohr wandern. 
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Abb. 25, Entwicklung der laminaren Geschwindigkeit in Rohrmitte an der vorderen Turbulenzfront. Re = 2550. 
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Abb. 26. Intermittenzfaktor als Funktion der Lauflange. Nach geglitteten Kurven von Abb. 19. 


9. Zusammenfassende Diskussion der Ergebnisse. Aus den mitgeteilten Untersuchungen geht 
hervor, daB in der Nahe der kritischen Reynoldszabl abwechselnd turbulente und laminare 
Partien bei zeitlich konstanter DurchfluBmenge durch das Rohr hintereinander herstrémen. 
Dabei bildet sich iiber den Querschnitt jeweils die der ausgebildeten turbulenten bzw. laminaren 
Strémung entsprechende Geschwindigkeitsverteilung aus. Die Verteilung hinsichtlich der Langen 
dieser Partien ist unregelmaBig. Bei der bisher im Schrifttum vertretenen Ansicht sind zeitliche 
Anderungen der DurchfluBmenge die entscheidende Ursache fiir das intermittierende Auftreten 
turbulenter Str6mung im Rohr’. Da bei den hier besprochenen Versuchen diese DurchfluB- 


1 Siehe FuGBnote 3 von S. 258 und Fuf8note 1 von S. 267. 
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mengenanderungen ausgeschlossen waren, ergibt sich, daf§ die eingangs erwahnte Erklarung 
nicht zutrifft. 

Uber das wirkliche Wesen dieser Strémung geben am besten die bei verschiedener Rohrlange 
versuchsmaBig bestimmten Werte des Intermittenzfaktors und der Intermittenzzahl n Auf- 
schlu8. Es handelt sich um die sukzessive Entwicklung voll ausgebildeter turbulenter Stromung. 
Dies wird deutlich, wenn man den Intermittenzfaktor y iiber der Lauflinge x fir konstante 
Werte der Reynoldszahl auftraigt, Abb. 26. Diese Darstellung wurde aus Abb. 19 abgeleitet, 
wobei fiir y die Werte der geglatteten Kurven genommen wurden. Der Intermittenzfaktor steigt 
stetig mit der Lauflange an, bis er den Wert 1 erreicht, d. h. bis voll turbulente Strémung eta- 
bliert ist. Dieser Vorgang erstreckt sich uber sehr groBe Rohrlangen, die fiir kleine Reynolds- 
zahlen sogar nach Tausenden von Durchmessern zahlen. 
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Abb, 27. Dimensionslose Intermittenzzahl als Funktion der Lauflainge. Nach geglatteten Kurven von Abb, 20. 


Abb. 27 zeigt den Verlauf der dimensionslosen Intermittenzzabl n d/U itiber der Rohrlange, 
der von den geglatteten Kurven der Abb. 20 umgezeichnet wurde. Die Intermittenzzahl nimmt 
im allgemeinen mit der Rohrlinge ab. Physikalisch bedeutet das Anwachsen von y und das 
gleichzeitige Abnehmen von n, daf die turbulenten Partien sich ausbreiten; sie wachsen und 
zehren die zwischen ihnen liegenden Jaminaren Partien sozusagen auf. Nach und nach schlieBen 
sich so die turbulenten Partien zusammen, so da dadurch die Intermittenzzah] abnimmt. 

Den unterschiedlichen Verlauf der Intermittenzzahl n tber der Rohrlange bei verschiedenen 
Reynoldszahlen kann man sich erklaren, wenn man sich aus den Werten des Intermittenzfaktors 
und der Intermittenzzahl die durchschnittlichen Vorbeistrémzeiten der Partien errechnet. Der 
Kehrwert der Intermittenzzah! n stellt die durchschnittliche Dauer t vom Beginn einer laminaren 
Partie bis zu Beginn der nachsten laminaren Partie dar. Auf Abb. 28 ist diese Dauer in dimen- 
sionsloser Form ¢ U/d iiber der Reynoldszahl aufgetragen, wie sie sich aus den Auszahlungen 
der Hitzdrahtschriebe, Abb. 11, ergibt. AuSerdem kann man auch die durchschnittliche Durch- 


gangszeit einer turbulenten Partie ermitteln zu 


gw 
deed os 


(34) 


Auch dieser Wert ist auf Abb. 28 mit eingezeichnet. Die Differenz zwischen den beiden Kurven 
entspricht der durchschnittlichen Dauer der laminaren Partien 


ei, (35) 


Die dimensionslose Dauer tU/d hat im vorliegenden Fall bei ungefahr Re = 2500 ein sehr 
ausgepragtes Minimum; fiir gréBere und kleinere Reynoldszahlen geht ¢U/d sehr schnell gegen 
unendlich. Die durchschnittliche Dauer der turbulenten Partien sind bei kleinen Reynoldszahlen 
klein und steigen mit Re sehr rasch an. Die durchschnittliche Dauer der laminaren Partien ist 
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bei kleinen Reynoldszahlen groB und nimmt ab, wenn Re wachst. Unter der Voraussetzung, 
dafi die Laufgeschwindigkeiten V der Turbulenzfronten etwa gleich der mittleren Durch- 
fluBgeschwindigkeit U sind, stellen die GréBen Tr U/d und t,U/d die durchschnittlichen 
Langen der turbulenten bzw. laminaren Partien als Vielfaches des Rohrdurchmessers dar. 


Da die Abnahme der Intermittenzzahl 
durch Zusammenschmelzen der turbulenten 
Partien entsteht, ist dieAnderung von n immer 
dann groB, wenn die laminaren Partien im 
Durchschnitt kurz sind. Das ist der Fall, wenn 
der Intermittenzfaktor y groB ist, also bei den 
groBeren Reynoldszahlen. Sind die laminaren 
Partien im Durchschnitt dagegen lang, wie es 
bei den kleineren Reynoldszahlen ist, so ver- 
schmelzen nur selten zwei turbulente Partien 
miteinander. Daher dndert sich n dann sehr 
wenig oder gar nicht, wie aus Abb. 27 hervor- 
geht. 

Aus dem Anwachsen von y und der Ab- 
nahme von n kann man quantitative Anhalts- 
punkte fiir die Ausbreitungsgeschwindigkeit 
der turbulenten Partien gewinnen. Ausbreitung 
einer turbulenten Partie hei®Bt, daB die Lauf- 
geschwindigkeit V, (bezogen auf das ruhende 
System) ihrer vorderen Front gréBer als die 
Laufgeschwindigkeit V, der hinteren Front 
ist. Passiert eine turbulente Partie von der 
Dauer try eine Stelle x des Rohres, so ist ihre 
Dauer an der Stelle x + dx 


1 iL 
ltr = === : 36 
tr a dtr tr 4 (7, A dx ( ) 
Da nun in der Zeiteinheit n Partien vorbei- 
stromen, so wird die Zunahme des Intermittenz- 
faktors 


eS ha 
dy =(5,— 7) 4- (37) 
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Abb. 28. Durchschnittliche Dauer der turbulenten und laminaren 
Partien und Gesamtdauer einer >vollen Periode x nach Hitzdraht- 
schrieben als Funktion der Reynoldszahl. Rohreinlauf Nr. 1. 
Lauflange x/d = 322. 
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Abb. 29. Zur Ermittlung der Ausbreitungsgeschwindigkeit der turbulenten Partien. 


Daraus folgt fiir den Zusammenhang zwischen Intermittenzfaktor und Lauflange x die Beziehung 


x 


1 rn 


(38) 
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i is y al eine Integrationskonstante dar. Es liegt nun die Vermutung nahe, daB 
Eee ene se Einlauf, die von der GréBenordnung der Hanlan oe, der 
n mag, die durchschnittlichen Laufgeschwindigkeiten V von der 

Wenn man noch annimmt, daf die statistischen Korrelationen 
hlassigbar sind, so lat sich der Geschwindigkeitsfaktor in (38) Wi 


laminaren Rohrstrémung sel 
Lauflinge « unabhangig sind. 
zwischen V,, V, und n vernac 


vor das Integralzeichen setzen: 
x 


il I. = P 
ayy chy Lee aie (39) 
gets ie 7) 


Auf Grund dieser Uberlegungen wurde auf Abb. 29 der Intermittenzfaktor tiber dem Integral | ; 


ete aufgetragen. Die Ergebnisse fiir konstante Reynoldszahlen lassen sich fiir x/d > 100 | 
U 


als gerade Linien darstellen. Die Laufgeschwindigkeiten der Turbulenzfronten scheinen danach 1 
tatsichlich von der Lauflange x unabhangig zu 


og sein 

a7 nan ; 2 

ri Aus der Neigung der Geraden kann man die 
an We relative Ausbreitungsgeschwindigkeit 


U U_ V,—Vv, U? 
Ve Cee ) 
ermitteln. Diese ist in Abb.30 als Funktion der |) 
Reynoldszahl dargestellt; sie darf wohl als eine ip 
Eigenschaft ausgebildeter Rohrstrémung be- 
trachtet werden, die von der speziellen Aus- 
006 bildung des Rohreinlaufs unabhangig ist. Nach 
Abb. 30 fallt die Ausbreitungsgeschwindigkeit 
mit abnehmender Reynoldszahl sehr schnell ab. |) 
006 Die Reynoldszahl, bei der sie den Wert Null er- | 
reicht, muf als die wirkliche kritische Reynolds- 
zahl angesehen werden, unterhalb derer sich die 
turbulente Rohrstrémung nicht mehr aufbauen 
kann. Der genaue Wert dieser kritischen Rey- 


(40) | 


010 


G04 


002 | noldszahl kann aus vorliegenden Versuchen nicht 
) ye . o 
eg bestimmt werden; er mag etwa bei Re = 2000 
ke liegen. 
Uplate ah Sai oS ee oa LD : 
Oe 2200 2400 2600 In dem untersuchten Bereich von Reynolds- 
Abb, 30. Relative Ausbreitungsgeschwindigkeit der turbulenten zahlen ist die relative Ausbreitungsgeschwindig- 


Partien als Funktion der Reynoldszahl. keit klein, d.h. die Differenz zwischen den Lauf- 
geschwindigkeiten der vorderen und _ hinteren 
Turbulenzfront ist klein gegen die DurchfluBgeschwindigkeit U. Die im vorigen Abschnitt auf 
Grund des VerhaJtens der Str6mung an den Turbulenzfronten ausgesprochene Vermutung, 
da} die Laufgeschwindigkeiten der Turbulenzfronten nicht viel von der DurchfluBgeschwin- 
digkeit abweichen, ist mit diesem Befund im Einklang. Zur Untersuchung der Laufgeschwin- 
digkeiten der Turbulenzfronten sind weitere Forschungsarbeiten notwendig. Diese wiirden 
nicht nur zur Klarstellung der Verhaltnisse bei der Rohrstrémung beitragen, sondern lassen 
dariiber hinaus eine wichtige Vertiefung der Kenntnisse iiber die grundlegenden Zusammen- 
hange bei der Ausbreitung turbulenter Stromungen erwarten. 

Aus den Untersuchungen ergibt sich nun das folgende Bild fiir die Etablierung turbulenter 
Rohrstrémung in der Nahe der kritischen Reynoldszahl: Durch die geometrischen Verhaltnisse 
am Einlauf (und natiirlich auch durch den Turbulenzgrad der zustrémenden Flissigkeit) werden 
im Rohr einzelne Partien turbulenter Strémung erzeugt, die nach Lange und Haufigkeit statistisch 
unregelmaBig verteilt sind. Auf dieses Phainomen, iiber das schon L. Schiller} einige Beob- 
achtungen mitgeteilt hat, soll hier nicht weiter eingegangen werden. Der Rohrstrémung ist es 
dann selbst tiberlassen, aus dieser Anfangsverteilung den voll turbulenten Zustand aufzubauen. 
Die Partien turbulenter Strémung werden mit der Strémung fortgefiihrt; dabei wachsen sie 


' Siehe FuBnote 1 von S. 258, 
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mit einer von der Reynoldszah] abhingigen Ausbreitungsgeschwindigkeit und schmelzen nach 
und nach zusammen, bis schlieBlich ein voll turbulenter Strémungszustand hergestellt ist. Diese 
Entwicklung kann sich iiber sehr groSe Rohrlingen erstrecken. Dieser Tatsache ist es zuzu- 
schreiben, da sich bei Messung zeitlicher Mittelwerte (Druckabfall, Geschwindigkeitsverteilung) 
der Wechsel von laminarer in turbulente Strémung stets als allmablicher Ubergang zeigt, der 
sich iiber ein endliches Intervall von Reynoldszablen erstreckt. Dies ist selbst dann der Fall, 
wenn grobste Mittel zur Erzeugung turbulenter Strémung angewandt werden. 

Es zeigt sich hier eine interessante Parallele zur Entstehung der turbulenten Grenzschicht 
an der Platte’. Bei dieser beginnt die turbulente Strémung — ebenfalls statistisch unregelmaBbig 
— an einzelnen Stellen. Die entstandenen ,,Turbulenzflecken® breiten sich aus, wobei sie mit 
der Str6mung fortgetragen werden, und iiberdecken schlieBlich die ganze Flache. Die Ausbreitung, 
die an der Platte zwei-dimensional vor sich geht, stellt sich im Rohr ein-dimensional dar. 

Eine praktische SchluBfolgerung, die man fiir den Begriff der Ubergangs-Reynoldszahl beim 
Rohr hieraus zu ziehen hat, ist, daB ein Zusammenhang zwischen Lauflange des Rohres und der 
Reynoldszahl, bei der die Strémung voll turbulent ist, besteht. Diese Erkenntnis hat schon 
Schiller niedergelegt. Wenn man genaue Aussagen machen will, unter welchen Umstinden in 
einem Rohr bei gegebenen Einlaufbedingungen turbulente Strémung zu erwarten ist, so geniigt 
nicht die Angabe einer auf den Rohrdurchmesser bezogenen Reynoldszahl, sondern es mu auch 
eine Rohrlange mit genannt werden. 

Nach den vorliegenden Erkenntnissen méchte man bei der Entstehung turbulenter Strémung, 
sei es im Rohr, in der Grenzschicht eines umstrémten Kérpers oder sonst wie, drei Phasen 
unterscheiden. Die erste Phase ist die Entstehung anwachsender Wellen, der Tollmienschen 
Wellen — infolge in der ankommenden Strémung vorhandener Stérungen, Schallwellen, Erschiitte- 
rungen oder dergleichen —, oder die Aushildung von WirbelstraBen — bei nicht stromlinien- 
formigen Kérpern —. In der zweiten Phase gehen diese Wellen oder diskreten Wirbel in 6rtlich 
begrenzten Gebieten in ungeordnete Stré6mung uber, wie man sie unter dem Begriff Turbulenz 
versteht. In der letzten Phase breiten sich dann die Gebiete turbulenter Strémung aus, bis sie 
den ganzen, in Frage kommenden Raum ausfiillen. Die vorliegenden Untersuchungen Jieferten 
fiir die Robrstr6mung einen Beitrag zu der letzten Phase der Turbulenzentstehung. Eine Anzahl 
von Fragen blieb dabei offen; iiber andere konnten nur Vermutungen ausgesprochen werden. 


Es ware eine lohnende Aufgabe, durch weitere Untersuchungen die verbliebenen Liicken zu 
schlieBen. 


(Eingegangen am 25. Januar 1956.) 
Anschrift des Verfassers: J. Rotta Gittingen, Bunsenstrafe 10. 


u . Aero. Sci. 18, (1951); S. 490. 
& é “a erdies a Ss. ne aa eee on the mechanics of boundary-layer transition, 
NACA TN 3489 (1955). 
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Zur Berechnung der kritischen Geschwindigkeit der reibungs- und 


relaxationsfreien thermochemischen Gasstrémung durch Rohr und Diise | 


Von H. Stiimke 


1. Einleitung. Bei der Auslegung von Strahltriebwerken interessieren u. a. folgende Fragen: 
Erstens, bei welcher Energiezufuhr bzw. bei welchem Querschnittsverhaltnis des Diffusors 
thermische Verstopfung der Brennkammer eintritt!»?, zweitens, bei welchem Querschnittsver- 


haltnis der Feuerdiise isentropische Entspannung auf Uberschall méglich ist®. Ganz ahnliche | 


Fragestellungen kommen auch in der Theorie des ebenen DetonationsstoBes vor**®. Im klassi- 


schen Fall eines Gases konstanter spezifischer Warme tritt der kritische Zustand gerade bei 6rt- | 


licher Schallgeschwindigkeit ein. Da dieser Satz unter gewissen Voraussetzungen auch fir die 
chemisch veranderliche Strémung gilt, ist fiir den Fall der Disenstrémung ohne Energiezufuhr 
wie fiir den Fall der Rohrstr6mung mit Energiezufuhr zwar bewiesen *°, jedoch scheint dies 
weniger bekannt zu sein. Wir wollen daher im folgenden die kritischen Geschwindigkeiten 
fiir diese beiden Falle noch auf einem etwas anderen Wege direkt berechnen, ihre Uberein- 
stimmung mit der 6értlichen Schallgeschwindigkeit zeigen und schlieBlich die numerische Be- 
rechnung an einem Beispiel erlautern. 


Unsere Annahmen sind: Reibungslosigkeit der Strémung, stationare Verhalinisse, Anwend- | 
barkeit der elementaren Stromfadentheorie und relaxationsfreie Einstellung der thermischen 


und chemischen Gleichgewichte. 

Auf Grund der letzten Voraussetzung ist der thermodynamische Zustand — ebenso wie bei 
konstanter chemischer Zusammensetzung — durch jeweils zwei GréBen voéllig bestimmt. Fir 
rein theoretische Betrachtungen eignen sich sehr gut der Druck und die Dichte als Zustands- 
variablen, da man dann die Temperatur als ZustandsgréBe ganz entbehren kann®. Da wir jedoch 
weiter unten eine spezielle kritische Geschwindigkeit numerisch berechnen wollen und die dafiir 
benétigten StoffgréBen im allgemeinen als Funktionen von Druck und Temperatur tabuliert 
vorliegen, wahlen wir die letztgenannten GroBen von Anfang an als Zustandsvariablen. Anderer- 
seits ist dieser Ansatz allgemein genug, um die im folgenden angewandte direkte Methode zur 
Berechnung kritischer Zustande bei beliebiger Wah] der Zustandsvariablen darzulegen. 


Haufig wiederkehrende Symbole sind: 
p = Druck; p als Index: A, = (5), Ausnahme: ¢,,C.. 6. tl. 
o = Dichte; 
T = Temperatur; T als Index: Arp = (=n) $ 
oT), 
cal,; 


2 
—— = 8314,4 one : 
MOlChem, Sk. K 


R = allgemeine Gaskonstante = 1,9865 
s? mol chem. sk. & 


M = Molekulargewicht; 
n = Konzentration im chemischen Gleichgewicht, als Molenbruch; 

Cp, h, s = spezifische Warme bei konstantem Druck, Gesamtenthalpie und Entropie pro Massen- 

einheit; 
C,, H, S = dieselben GréBen pro Mol; 

S* = Entropie pro Mol beim Druck von 1 Atm.; 
C = Molwarme des Gasgemisches bei konstantem Druck und konstanter Konzentration; 
F = Rohr- bzw. Diisenquerschnitt; 
w = Strémungsgeschwindigkeit ; 


1 E. Sanger, Madras Institute of Technology 3 (1954) S. 91. 
2 I. Sanger-Bredt, Bemessungstafeln fiir Staustrahlantriebe. Bericht Nr. 4, Teil 1, des Forschungs- 
institutes fiir Physik der Strahlantriebe e.V., Stuttgart Gm Druck). 

5 E. Sanger, Schweizer Arch. angew. Wissensch. u. Technik, 16 (1950) S. 341 und S. 369. 

4 K. Oswatitsch, Gasdynamik, Wien 1952. 

5 R. Becker, Z. Physik 8 (1922) S. 321. ; 


6 J.O. Hirschfelder, C. F. Curtiss und R. B. Bird, Molecular theory of gases and liquids, New York 1954. 
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Wr = kritische Geschwindigkeit der Strémung im Rohr; 

wy = kritische Geschwindigkeit der Strémung in der Lavaldiise; 
a = Schallgeschwindigkeit ; 
i = Index fiir eine chemische Komponente des Gasgemisches ; 
0 = Index fiir Verhiltnisse im Anfangsquerschnitt; 


2. Schallgeschwindigkeit. Die Scha'lgeschwindigkeit ist im allgemeinen frequenzabhangig. 
Wir betrachten zunachst adiabatische Zustandsanderungen, die so langsam erfolgen, daB man 
die Einstellung der thermischen und chemischen Gleichgewichte als verzigerungsfrei ansehen 
kann. Die Zustandsinderung erfolgt dann isentropisch: 

s(T, p) = konst. (1) 
und ihre Fortpflanzungsgeschwindigkeit (Schallgeschwindigkeit a) ist durch die Formel 


gegeben. In dieser Gleichung kénnen wir auf Grund unserer obigen Annahmen die Dichte o als 
Funktion von p und T allein ausdriicken: 


@ =o(p; T). (3) 
Ein zeitlich infinitesimal benachbarter Zustand eines Masseteilchens unterscheide sich vom Zu- 


stand T’, p durch die Variationen dT, dp, do. Fir diese gelten dann die beiden Differentialbe- 
ziehungen 


spdp + spdT = 0 (1) 
und 
Op dp + er dT —do = 0. (3)’ 
Aus ihnen folgt unmittelbar fiir die Fortpflanzungsgeschwindigkeit a sehr langsamer Zustands- 
anderungen die Formel 
op ay 12 Se 
2 = = = 3 4 
ys (3), Op ST — OT Sp @ N, (4) 


wobei zur Abkiirzung 
aI’ 
Na =) (ep ST — QT Sp) 
gesetzt wird. 

Die Gesamtenthalpie h des Gasgemischs mége nun so definiert sein, daB sie auBer dem 
Produkt pv noch die gesamte innere Energie einschlieBlich der latenten chemischen Energie 
enthalt. Nach dem ersten und zweiten Hauptsatz der Warmelehre ist dann bei reversibler 
W armezufuhr 


Ure, = T ds = dh ——- dp. (5) 


Fa8t man auch s und h als Funktionen von p und T auf, so folgen aus (5) sofort die beiden Be- 
ziehungen 

TE sT = hr Py (6) 

ae he es: (7) 

(& 
Durch Einsetzen derselben in (4) erhalt man fiir die Schallgeschwindigkeit a sehr langsamer 
Zustandsanderungen noch die Darstellung: 
hy 


ee (8) 


2 reece 
ei ep hr —or(hp —-] 


Bei schnelleren Zustandsanderungen, also z. B. bei gewéhnlichen Schallwellen mit héheren 
Frequenzen, kann der Fall eintreten, daB die Einstellung der thermischen Gleichgewichte noch 
relaxationsfrei erfolgt, wahrend die trageren chemischen Gleichgewichte nicht mehr realisiert 
werden; die chemische Zusammensetzung bleibt dann beim Durchgang der Schallwellen prak- 
tisch konstant. In einem solchen Frequenzbereich kann man die obigen Formeln (4) bzw. (8) 
unter Festhalten der Konzentrationen anwenden. 
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Es ist auch der Fall denkbar, da die Einstellung der Dissoziationsgleichgewichte nur in 
Anwesenheit eines Katalysators vonstatten geht, und da man diesen nach Belieben hinzufiigen 
oder entfernen kann. In Abwesenheit des Katalysators miissen dann sogar bei sehr langsamen 
Zustandsinderungen die Konzentrationen in (4) bzw. (8) als konstant behandelt werden. 


3. Kritische Geschwindigkeit in der Lavaldiise. In der Lavaldiise sind bei Abwesenheit von Warme- 
zufuhr, Reibung und StoBverlusten nicht nur die Zustandsgleichung (3) und die Entropieglei- 
chung (1) gewahrt, wobei sich jetzt die raumliche Koordinate x in Stromrichtung als natiirlicher 
Parameter anbietet, sondern es gilt auch noch die Energiegleichung 


z + A(T, p) = konst. (9) 
sowie die Kontinuitatsgleichung 
F 9 w = konst. (10) 
Im engsten Querschnitt der Diise erreicht die Stromdichte ein Maximum: 
0 w = max fiir F = Finin« (11) 


Im allgemeinen nehmen dort auch alle anderen ZustandsgréBen einen Extremwert an. Sofern 

dies nicht der Fall ist, pflegt man von einem ,,kritischen‘‘ Zustand zu sprechen!. Auf jeden Fall 

werden in einem Diisenquerschnitt, welcher der engsten Stelle infinitesimal benachbart ist, fur 

die Variationen dp, dT, dg und dw aufer den Gleichungen (1)’ und (3)’ noch die Beziehungen 

gelten 

wdw-+h,dp +hrdT =0, (9)’ 

odw+wdo=0. (1) 

Die Gleichungen (1)’, (3)’, (9)’, (11)’ bilden ein System von vier linearen homogenen Glei- 

chungen. Beim Vorliegen eines Extremfalles werden die Gleichungen durch identisch verschwin- 
dende Variationen erfiillt. Im kritischen Zustand dagegen muf die Determinante 


w hy hr 0 
0 0 0 w 


OSe> ser 1 
Onsen sree) 


=0 (12a) 


sein. Subtraktion der mit T multiplizierten vierten Zeile von der ersten Zeile und Addition der 
mit w multiplizierten dritten Zeile zur zweiten Zeile fiihrt unter Beachtung der Gleichungen (6) 
und (7) zu der vereinfachten Form 


fw oe 0 
; 0 12b 
Q wep worl” a 
0 Sp ST 
aus der sofort 
wy = ——T . (12c) 


Op ST — OT Sp 
folgt. In Verbindung mit (4) ergibt sich hiernach, da® unter den einleitend genannten Voraus- 


setzungen die kritische Diisengeschwindigkeit wz, auch bei thermochemischer Gasstrémung mit 
der Schallgeschwindigkeit fiir geringe Frequenzen tibereinstimmt. 


4. Kritische Geschwindigkeit im Rohr. Beim Durchstrémen der Strecke x4, Xp eines Rohres 
werde die pro Zeiteinheit eintretende Masse auf das (1 + yw)-fache, der eintretende Gesamt- 
impuls auf das (1 + »)-fache und die eintretende Gesamtenergie auf das (1 + e)-fache vergréBert. 
Dabei wollen wir uu, 7 und ¢ als unabhangig vom Endzustand in x voraussetzen. In der iiblichen 
Behandlung der Rohrstrémung wird meistens 4 =n =€ = 0 gesetzt. Am Anfang der Strecke 
oder nach erfolgter Massenzufuhr wird i. a. ein chemischer Nichtgleichgewichtszustand bestehen. 
Auf jeden Fall soll aber am Ende der Strecke thermochemisches Gleichgewicht erreicht sein, 


1 Siehe FuBnote 4, S, 282. 
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Man hat dann fiir die Zustande in Punkten X) <«x4und x > xp auBer der Zustandsgleichung (3) 
die Beziehungen 


52 
_ Energiegleichung: - eC fh) & 4 he] ere), (13) 
Impulsgleichung: P + 0 w = (po + Oy w?) (1 +n), (14) 
Kontinuitatsgleichung: ow=ew(1+ 4). (15) 


Die rechts stehenden GréBen py, 0) usw. hangen im allgemeinen von mehreren Parametern 
ab. Wir nehmen jetzt fiir alle diese Parameter aufer einem, der y hei®en mége, konstante Werte 
an. Die GréBen po, @) usw. sind dann Funktionen von y. 


Eliminiert man jetzt mittels der Zustandsgleichung (3) und der Kontinuitatsgleichung (15) 
die GréBen w und ¢ aus der Energiegleichung (13) und der Impulsgleichung (14), so stellen die 
Energiegleichung (13) und die Impulsgleichung (14) in der p, T-Ebene je eine Schar von Kurven 
f(p, T, y) = 0 bzw. g(p, T, y) = 0 mit stetig sich drehender Tangente dar. Reelle Schnittpunkte 
zweier Kurven f und g mit gleichem Parameter y bedeuten mégliche Zustande im Punkt x. 
Der Grenzfall gegen die physikalisch nicht realisierbaren Zustande tritt offenbar dann ein, wenn 
sich die Kurven f und g beriihren. Formuliert man diese Bedingung analytisch, so erhalt man eine 
Beziehung zwischen der kritischen Geschwindigkeit wr der Rohrstrémung und den Zustands- 
variablen p und T. Dies ist das tbliche Verfahren, mit welchem man auch bei Annahme eines 
(p, @)-Systems die kritischen Zustande zu bestimmen pflegt!. 

an kann nun, ohne den mathematischen Gehalt der Herleitung zu andern, den Beweisgang 
wesentlich abkiirzen, indem man direkt bei konstantem y eine nichttriviale Lésung fiir die Dif- 
ferentiale der ZustandsgréBen sucht. 

Zu den friiheren Gleichungen (3)’, (9)’, (11)’ tritt dann noch die aus der Impulsgleichung 
folgende Relation 


20wdw-+dp+w'doe =0. (14)’ 


Als Parameter, nach dem die Differentiation ausgefiihrt wird, wird man bei dieser Deutung eine 
der Zustandsvariablen w, 9, p oder T nehmen. Eine nichttrivale Lésung des Problems kann 
nur dann bestehen, wenn die Determinante 


w ih 
Zio) «1 Ors 40" 
OPT Op te al 
0 Q —W. @ 


= 0 (16 a) 


ist. Subtraktion der mit w multiplizierten vierten Zeile von der zweiten Zeile und Addition der 
mit w multiplizierten dritten Zeile zur vierten Zeile fihrt zu der einfacheren Form 


w hp hr 
ow 1 0 |=0, (16b) 
Q@ WOp WOT 


aus der sich sofort 


tee ba (16c) 


ergibt. Ein Vergleich mit (8) zeigt, daB unter den oben genannten Voraussetzungen die kritische 
Geschwindigkeit wr des Rohres auch bei dissoziierenden Gasen mit der Fortpflanzungsge- 
schwindigkeit a niederfrequenter Stérungen tibereinstimmt. 


5. Spezielle thermochemische Ansatze. Fiir die Durchrechnung eines Zahlenbeispiels be- 


| nétigen wir einen ausfiihrlichen Ansatz fiir die Zustandsgleichung eines reagierenden Gasge- 


misches. Da hieriiber noch wenig bekannt ist, setzen wir voraus: Die Temperatur sei so hoch 


» bzw. der Druck so gering, daB die Zustandsgleichung der Einzelgase und des Gasgemischs in der 


1 Siehe FuGnoten 4 bis 6 auf S. 282 
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Form ee 
aca 17 
Chia Re (17) 
angenommen werden kann. 
Dabei wird ar 
60 p Mp Bee ke Tr), 17)! 
eo == (1 + Pa), or =p (1 M (17) 
Zwischen den Funktionen h, s einerseits und H, S andererseits bestehen dann die Beziehungen 
H S 
= ae 18 
i ive S=ar° (18) 
Hier ist! 
M = 5 Myn,, n= Dns (19) 


C= Gt H= 3) Hi; n;, S= Sin, (20) 


die 
H, = [ C,:dT +k, 


C dT C dT eh) 
= |B peediod — . . = Pi — — . . 
8 =[ ip Ring +h = | T Rln p— Rlnn; +1;. 


i Ty 


Fiir die spitere Rechnung werden wir noch die partiellen Ableitungen der Funktionen H 
und S nach T und p brauchen. Sie lauten in dimensionsloser Schreibweise 


Ar cae (G 7 H; T on; TS en ‘S; T on; 
R R Sipe luo Dee ie AR af Ta sO ars 
u L (22) 
pAp _ H; pon; P Sp shee Si p On; 
peed Ge Py wae TREN Aye TR tals DR op 


Die Konstanten k; in (21) mégen irgendwie so normiert sein, dafs die Funktion h die Bedeu- 
tung der Gesamtenthalpie im Sinne des Abschnitts 2 annimmt. Entsprechendes soll fiir die 
Konstanten /; in der Definition der Entropie gelten. 


6. Zahlenbeispiel: Kritische Geschwindigkeit von dissoziierendem CO,. Zur Erlauterung der 
Rechnung betrachten wir die Verhaltnisse, die bei der Strémung von dissoziierendem Kohlen- 
dioxyd beim Druck von einer physikalischen Atmosphare und einer Temperatur von 2000° Kauf- 
treten. Als Reaktionsgleichung nehmen wir an: 


2CO, =2CO0+0,. (23) 


Bezeichnet « den Dissoziationsgrad des Kohlendioxyds im Gleichgewicht, so sind die Molen- 
briiche 


of 
l—«a o a 
nco, = Ao Pol Oa we (24) 
1+— 1+— 1++ 


Nach Ulich? besteht zwischen diesen die aus dem Massenwirkungsgesetz folgende Beziehung: 


a3 
PCo Po, Néo No, 2) 
= pe ae SNe a : (25) | 
Co, CO, (: ue rae bs 


Dabei ist die Funktion K, als reine Temperaturfunktion mit den Werten von Tabelle 1 gegeben, | 
zu denen wir gleich die Differenzen bis zur dritten Ordnung berechnen 


1 I. Sanger-Bredt, Z. angew. Math. Phys. 6 (1955) S. 35. 
* H. Ulich, Kurzes Lehrbuch der physikalischen Chemie 4. Aufl. S. 111/2, Dresden und Leipzig 1942. 
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Tabelle 1. Daten zur Berechnung von o und Li log aN 
oT Atm 
fiir p = 1 Atm., T = 2000 °K 
K 
T [°K] log ic A A As 
wor | stat | on 
Sa nag Se A ee a 2,95 
2.000 — 5,75 2.95 — 1,92 0.83 
2500 — 2,80 1,86 — 1,09 
3000 — 0,94 


Gleichung (25) stellt implizit « und damit die Konzentrationen n; als Funktionen von p und T 
dar. Durch partielle Differentiation von (25) erhalt man die Relationen 


a 
zx a(1 eg | ee 


Clea 
dp Pp 3 ; 
25) 
2(1 +z)a—9 (5) 
Ou 2 @ Noo Ke 
aT 3 aT \'°8 Atm 
Durch entsprechende Differentiation der Gleichungen (24) folgt 
3 3 
Onco, | Oa 2 Anco, __ 0a 2 
oT oT (: tel dp Op (: oe 
2 2 
Anco ox 1 Ono Ox 1 
2 2 
i 
Ono, 0x Eos Ong, Ae 2 
oT oT ( +3) dp Op (: ei) 
Zi 2 
_- Zusammenfassung der Gleichungen (24)’ und (25)’ gibt 
TOnco, —a(1—a) Ta (log Ke Pn, a (1—a) 
oT 9 ( =) oT Atm. Op 21 =) 
2 2 
T On¢o a(l1—a) To Kp P ONO —a(]—«) 
Pinar aon? Rea” fe = ae (26) 
3 ( et a 2 3(1 + zs 
2 2 
oo, eV) oa =) pam,  —a(l—a). 
oT 6(1 oa oT Atm. Op 6(1 = 
7 ” 
Fir p = 1 Atm., T = 2000 °K ist nach Tabelle 1 log — = — 5,75. Durch Iteration errechnet 
sich daraus mittels (25) i 
ao = 0,015 1478. (27a) 


Ferner ergibt sich aus dem Differenzenschema der Tabelle 1 unter Benutzung der Stirlingschen 
Formel fiir numerische Differentiation: 


{l 1 1 : 1 
F (——)+ 6F (— a 1 OF 
Gira 7) +9 Gz) 1 ( 7) (3) 


du 2 6 2 


—- +: 
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Tabelle 2. Berechnung der Gemisch-Funktionen fiir p = 1 Atm., T = 2000 °K 
Mischgas Einzelgase 
GréBe Zahlenwert co, co 0, GréBe 
bzw. Summe 
1 44,010 28,010 32,000 M; [g/mol 
Ge 
2 7,2580 4,3590 4,5436 oe 
i; 
3 6,0620 20,7161 4,0820 aa 
S* 
4 37,1860 31,1030 32,3161 a 
5 Zn; 1 0,9774490 0,0150340 0,0075170 nj 
6 0,0228 4,1974 4,8906 ing 
S: 
a 37,2088 35,3004. 37,2067 R 
8 M [g/mol] 43,6791 43,0175 0,4211 0,2405 M; n; [g/mol] 
C CE 
9 7 7,19400 7,09432 0,06553 0,03415 oe 
H Ai; 
10 ae 6,26742 5,92530 0,31144 0,03068 nae 
s* S* 
11 R 37,0579 36,3474 0,4676 0,2429 RN 
12 = —n,ln n; 0,1222 0,0223 0,0631 0,0368 —n, In n; 
S S; 
13 R 37,1801 36,3697 0,5307 0,2797 Fee 
T on; oP On; 
14 aT 0 —0,10645661 | 0,07097107 0,03548554 aT 
15| TMrp[g/mol] |—1,56172 —4,68516 1,98790 1,13554 M; rae [g/mol] 
H; Ton; . 
16 sue Ge. 0,969755 | —0,645340 1,470243 0,144852 a ae 
ih T on; * 
17 ap | |-0,60454 —3,95870 2,20741 1,14675 “ es 
at T On; T On; 
cee | = n 
18 Za nn; ar | 0,46902 0,00242 0,29789 0,17355 Shiga 
S; T On; . . 
19 RE —0,13551 —3,96112 2,50531 1,32030 2 at 
On: 
20 P , ao 0 0,007403102 | —0,004935401 | —0,002467701 P On 
Pp dp 
On; 
21)  p Mp [g/mol] 0,1086035 | 0,3258105 | —0,1382406 | —0,0789664 | m2 ; % [g/mol] 
7 p 
P 
22 aa —0,0674381 0,0448776 | —0,1022425 | —0,0100732 H; p On; 
RETO p 
S¥ p On; 
z8 aS R Op 0,0420395 0,2752918 | —0,1535058 | —0,0797465 S# p dn; 
t R_ dp 
p On; 
24 > on ni —0,0326156 | 0,0001688 | —0,0207159 | —0,0120685 =n poe 
i Pp 
Si p On; S; pon; 
25 R dp 0,0094239 0,2754605 | —0,1742216 | —0,0918150 = ue ‘ 
i Pp 
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Tabelle 3. Berechnung einiger Kontrollgréfen fiir p = 1 Atm., T = 2000 °K 


GréBe Zahlenwert 
iG H; Tan 
26 = | i i 
oa ag a 8,16376 
TeSys 6 S; Tan 
27 = : ae 05849 
Rin gfe pe R oT uy 
P Sp S; P On; 
28 ==—14 — — 0,9905761 
T Mr 
= 1544 
29 Ni 0,03575 
p Mp 
30 0,002486395 
M 
31 a 0,01558328 
RY MM 
S pM 
32 — - 0,092 444 42 
R M 
H TM, 
== — 0,2240878 
ve RT M 
34 eee 2 1,3293522 
R M 
35 MT. TSy S TMr 8.38784 
R R R M 
Sa ae a ge ES 8,38785 
R R RT M 
37 SD Cs Sea, — 1,0830205 
raw R R M 
38 Wap \ be eM. oy, — 1,0830214 
Bia! yoo RT RT M 


Tabelle 4. Berechnung der Schallgeschwindigheit fiir p = 1 Atm., T = 2000 AK 


| GroéBe 


S| RK. ( Reo ROW Con lors RM M 
MT MT. 
40 3 a? => es sr} (CR ) 
C C 
Al | : eed) 
42 pees RT |m? s~? 
0 A 
43 a ea | 
44, a {m s~*] 


EEE 


Zahlenwert 


7,28695 


TOMO TEC 


1,161446 


380703,8 
438219,4 
661,981 
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die Ableitung von log ae an der Stelle p = 1 Atm., T = 2000 °K zu 


d Kp eet, eee ee 27h 
ar (leg a = 0,007 19000°K+, (log SP) = 14, 380 00. (27b) 


Wir haben die GréBen « und - log ea mit 6 Ziffern angegeben, um den voraussichtlich schwa- 
chen Einflu8 der Dissoziation auf die kritische Geschwindigkeit mit Sicherheit herauszuholen, 
obwohl natiirlich aufgrund der Zahlenangaben in Tabelle 1 nur drei Ziffern gerechtfertigt waren. 


Die einzelnen Schritte bei der Berechnung der Geschwindigkeit sind in den Tabellen 2 bis 4 — 


festgehalten. Die thermodynamischen Funktionen der Einzelgase in Tabelle2 wurden den Tabellen 
von Wooley 1 entnommen. Dabei konnten alle Entropiewerte ohne Umrechnung verwendet werden, 
da sie bereits so normiert sind, daB sie am Nullpunkt der absoluten Temperatur verschwinden. 
Die Nullpunktsenthalpie von CO, und O, wurde gleich Null gesetzt, diejenige von CO mit 
66677 cal,;/molcn. sk. angenommen, um innerhalb der Rechengenauigkeit die Gleichungen (6) 


und (7) exakt zu befriedigen. Dieser Zahlenwert stimmt bis auf etwa 1°/)) mit Angaben der | 


neueren Literatur fiir die Nullpunktsenergien der Reaktionspartner von Gleichung (23) tber- 
ein?. In Tabelle 3 wird eine Kontrollrechnung fiir die Erfiillung der Gleichungen (6) und (7) ge- 
geben, siehe Zeile 35 bis 38. 

Nach Zeile 40 bzw. 41 der Tabelle 4 liegt das Quadrat der Schallgeschwindigkeit im Falle 


konstanter Konzentrationen um etwa 1°%, héher als im Falle veranderlicher Konzentrationen. | 
Man kann sich dies dadurch plausibel machen, daB die spezifischen Warmen im ersten Falle | 
niedriger sind und daher ihr Verhialtnis c, :c,, welches die Schallgeschwindigkeit maBgeblich | 


bestimmt, gréBer als im Falle veranderlicher Konzentrationen ist. Schallwellen, bei denen die 
chemischen Konzentrationen dauernd konstant bleiben, kénnen sich demnach in einer Diise 
entgegen der Strémungsrichtung durch die Zone der kritischen Geschwindigkeit hindurch fort- 
pflanzen. 

Fir das betrachtete Zahlenbeispiel ist ein geringer Dissoziationsgrad mit schwacher Tempera- 
tur- und Druckabhangigkeit der Konzentrationen charakteristisch. Um Gegenbeispiele zu bilden, 
wurde der Rechnung das Dissoziationsgleichgewicht von molekularem und atomarem Wasser- 
stoff bei 1 Atm. und 3500°K bzw. 5000°K zugrundegelegt. Im ersten Falle ist ein mittlerer Disso- 
ziationsgrad mit einer sehr starken Temperatur- und Druckabhangigkeit der Konzentrationen 
verbunden, im zweiten Fall ist ein sehr hoher Dissoziationsgrad wieder mit einer geringen Tem- 
peratur- und Druckabhangigkeit der Konzentrationen gekoppelt. Auch in diesen beiden Fallen 
errechnet sich fiir die Schallgeschwindigkeit bei konstanten Konzentrationen ein gréferer Wert 
als im Falle des Dissoziationsgleichgewichts. Die Abweichung betragt in diesen Fallen gréfen- 
ordnungsmafig 5%, also rund zehnmal so viel wie im obigen Beispiel. 


7. Zusammenfassung. Die kritische Geschwindigkeit der Rohr- bzw. Diisenstrémung eines 
chemisch veranderlichen Gasgemisches wird unter den Annahmen des permanenten chemischen 


Gleichgewichts, der Reibungslosigkeit und stationaérer Verhaltnisse nach der Stromfadentheorie | 


berechnet und mit der Schallgeschwindigkeit verglichen. Es zeigt sich, daB die kritischen Ge- 
schwindigkeiten mit der Fortpflanzungsgeschwindigkeit niederfrequenter Schallwellen iiber- 
eee in denen das chemische Gleichgewicht der Zustandsdnderung ohne Verzogerung 
olgt. 
Die Rechnung wird am Beispiel der Strémung von CO, bei 2000° K und 1 Atm. unter Dis- 
soziation in.CO und QO, erlautert. Hierbei ist die Schallgeschwindigkeit bei konstanter chemi- 
scher Zusammensetzung um etwa 0,5% héher als die kritische Strémungsgeschwindigkeit, 


was man mit dem héheren Wert der spezifischen Warmen im Falle des Dissoziationsgleichge- 
wichts erklaren kann. 


(Eingegangen am 26. Januar 1956.) 
Anschrift des Verfassers: Dr. H. Stiimke, Stuttgart-Degerloch, MeistersingerstraBe 27 


1 H.W. Wooley, The NBS-NACA tables of thermal properties of gases. Table 9.10. Molecular oxygen 


fone 3 ye Juli 1949; Table 14.10. Carbon monoxide (ideal gas state), Sept. 1951; J. Research NBS 52 


* H. Zeise, Thermodynamik III, 1, Tabellen, Leipzig 1954, 


‘Tragwerke aus Aluminium | 


Von Dr. sc. techn., LL. D. h. c., Dr.-Ing. h. c., Dr.-Ing. E. h. F. Stiissi, o. Professor an der Eidgen. 
Technischen Hochschule, Ziirich. Mit 174 Abbildungen. VII, 198 Seiten Gr,-8°. 1955. Ganzleinen 
DM 22,50 


Inhaltsiibersicht: I, Allgemeine Uberlegungen. — II. Der Baustoff und seine Eigenschaften: Die Her- 
stellung von Aluminium. Die Aluminiumlegierungen, Lieferformen. Festigkeit, Verformung und zulissige 
Beanspruchung. — III. Verbindungsmittel: Allgemeines, Nieten. Schrauben. SchweiBen. Kleben. — 
IV. Besondere Festigkeitsprobleme und Stabilititsprobleme: Biegung und Verdrehung. Knicken. Torsions- 
knicken und Kippen, Ausbeulen. — V. Ausbildung und Bemessung der Bauelemente: Zusammengesetzte 
Vollwandtrager. Fachwerktrager. Besondere Bauformen. — VI. Die Herstellung von Aluminiumtrag- 

_werken: Werkstattarbeiten, Montagearbeiten. Korrosionsschutz und Oberflaichenbehandlung. — VII. Aus- 
fiihrungsbeispiele: Briickenbauten. Hochbauten. Verschiedene Anwendungsformen. — VII. Obersicht 
iiber die Normen und mechanischen Eigenschaften der Knetlegierungen. 


_ Schon heute sind eine Reihe von bemerkenswerten Bauwerken aus Aluminiumlegierungen ausgefihrt | 
worden; die Bedeutung dieser Bauweise wird mit dem Ausbau der Produktionsstétten von Aluminium 
weiter ansteigen. Es scheint deshalb der Versuch erwiinscht, die wichtigsten Grundsitze fiir den Entwurf 
und die Ausbildung solcher Tragwerke aus den besonderen Baustoffeigenschaften abzuleiten. Wenn sich 
auch diese Grundsitze weitgehend an die des Stahlbaues anlehnen, so wire es doch verfehlt, die Bauformen 
und Bemessungsregeln, die sich im Stahlbap bewahrt haben, direkt. auf den Leichtmetallbau zu iiber- 

_tragen; es muB vielmehr angestrebt werden, die Bauformen hier so zu entwickeln, daB sie den Besonder- 
heiten des Baustoffes méglichst gut entsprechen. Dies bedingt vor allem eine’ geniigende Kenntnis von 
Herstellung und Eigenschaften des Baustoffes, wobei insbesondere die Frage-der Dauerfestigkeit eine ein- 
gehende Beriicksichtigung erfordert. Der niedrige Wert des Elastizitatsmoduls laBt gegeniiber Stahl die 
Bedeutung der Stabilitatsprobleme stark ansteigen, deren Darstellung deshalb ein verh4ltsnismaBig breiter 
Raum gewidmet ist. Ein Uberblick ber die wichtigsten bisher ausgefiihrten Anwendungsbeispiele soll die 
heute schon bestehenden Méglichkeiten der Aluminiumlegierungen im Bauwesen veranschaulichen. 


Statik der rahmenartigen Tragwerke 


Von Prof. Dr.-Ing. J. Pirlet, Kéln, ehem. Honorarprofessor der Technischen Hochschule Aachen, 


‘ Mit 80 Abbildungen und fiinf Tafeln in einer Tasche. VII, 168 Seiten Gr.-8°. 1951. 
Ganzleinen DM 24,— 


Inhaltsiibersicht: Erster Teil: Grundlagen der Berechnung rahmenartiger Tragwerke. ‘Der elastisch 
eingespannte Stab. Das Stabgefiige. Der Knoten. Der Knotenstabzug. Die Einspanngrade ¢. Berechnung 
der Einspannmomente M; und M,, d. h. der Unbekannten in einem Stabzugfelde l;. Gleichungen fiir M; 
und My, bei einzelnen Belastungsarten. = Anhang: Zusammenstellungen der Ergebnisse: — Zweiter 
Teil: Der elastisch eingespannte kontinuierliche Balken in rahmenartigen Tragwerken. — Das Berechnungs- 
yerfahren. Einleitung. Die Einspanngrade ¢ . Die Einspannmomente M; und My, des beiderseits elastisch 
eingespannten Balkens. Die. Kinspannmomente M in kontinuierlichen Trigern. — Dritter Teil. Zahlen- 
beispiel. Von der Belastung unabhingige Werte. Von der Belastung abhingige Werte. Vereinfachtes 


Verfahren. — Literaturverzeichnis. — Sachverzeichnis.. 
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EinfluBflichen fiir Kreazwerke 


plese . 
«Freiaufliegende und iiber mehrere Offnungen durchlaufende Systeme, Von Dr.-Ing. Hellmut Homberg, 
Hagen/Westf., und Dr.-Ing. Jos:f Weinmeister, Linz a./D. Zweite verbesserte Auflage. Mit 40 Ab- 


bildungen, VII, 156 Seiten 4°, 1956. : Ganzleinen DM 43,50 


Tahalteahereiekt: A. Erléuterungen zum Aufbau und Gebrauch der Lésungen. — Grundsitzliches — 
gar Ausbildung und Wirkungsweise von Kreuzwerken und verwandten Systemen. Begriffsbestimmung 
und Einfithrung. Bezeichnungen am Kreuzwerk. Beschreibung des Gleichungsaufbaus der Finflufflachen. . 
Hilfssystem in der Querrichtung. Lésungen und Tafeln fir die GréSen B;, baw. Bj,/,). Durchlaufende el a 
Kreuzwerke. Berechnang von Kreuzwaben (orthotropen Platten). Auswertung der EinfluBflachen. Hiil- | 
ites der Biegemomente und Querkrifte. Durchfiihrung der Berechnung, Naherungsweise Berechnung sc 
yon Kreuzwerken. Die Wirkung und naherungsweise Beriicksichtigung der Drehsteifigkeit. Die Wirkung BE 
der Schubsteifigkeit. Beispiele. — B. Lésungen fiir die EinfluBflichen von Kreuzwerken. Beidseitig frei = ; 
aufliegende Kreuzwerke iiber einer Offnung. Kreuzwerke iiler einer Offnung, links frei aufliegend, rechts 4 
fest eingespannt. Kreuzwerke mit unendlich vielen, unendlich schmalen Quertragern beiderseits starr 
eingespannt. Durchlaufende Kreuzwerke tiber zwei Offnungen. Durchlaufende Kreuzwerke iiber drei 
Offaungen. Ordinaten fry der Biegelinien des Balkens auf’zwei Stiitzen.:— C. Tafeln der Einflaiflachen 
usw. von Kreuzwaben (orthotropen Platten). — D. Allgemeine Lésungen fiir die ¢ Auflagerkriifte Bj, fie Balken 


auf drei bis zehn elastischen Stiitzen, — E. Tafeln der Auflagerkrafte Bix von Balken auf elastischen Stiitzen 


(Querverteilungszahlen). Balken auf 3 bis 8 Stiitzen. “Balken auf unendlich vielen Stiitzen. 


Es werden gebrauchsfertige, im Sinne der Baustatik genaue Lésungen fiir die EinfluBflachen der statischen 
GréBen geboten, die dazu dienen, dem Statiker eine schnelle Berechnung der mehrfach oder hochgradig 
statisch unbestimmten Kreuzwerke zu erméglichen, Die Neuauflage enthilt Lésungen fir freiaufliegende 

und iiber mehrere Felder durchlaufende Kreuzwerke und daher auch den im Vorwort zur foe Auflage © Me 
angekiindigten zweiten Teil des Gesamtwerkes. Es hat sich als zweckmaBig erwiesen, fiir durchlaufende 
Kreuzwerke solche Lésungen anzugeben, die beliebigen Stiitzweitenverhiltnissen entsprechen. Weiter 
wird ein Verfahren zur naherungsweisen Berechnung von Kreuzwerken angegeben. Die fiir die praktische 
Durchfithrung von Kreuzwerkberechnungen so niitzlichen Tafeln der Querverteilungszahlen wurden durch 
Aufnahme neuer Tafeln fiir die Randsteifigkeit 1 <r < 2 und fir den Balken auf oo vielen Stiitzen erheblick 
erweitert, sie kOnnen auch fiir das Naherungsverfahren verwendet werden, Die zugehérige Theorie ist inner- 
halb der Arbeiten »Kreuzwerke, Statik der Trigerroste und Platten‘ , Berlin 1951, und »-Beitrag zur 


Krelcworkiersekmaees Stahlbau 1954, verdffentlicht worden. 
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